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ABSTRACT 


Accurate enthalpies of fluid mixtures over wide range of pressure 
and temperature are directly useful in design and provide an improved 
basis for the theories of fluids, improvement of methods of prediction, 
and for derivation of other thermodynamic properties. The purposes of 
this work are: 
1) To construct a calorimeter for the existing recycle system to 
allow measurements on the isothermal effect of pressure on enthalpy. 
2) To obtain accurate data on the effect of pressure on enthalpy 
for mixtures of methane with carbon dioxide, with carbon dioxide and 
nitrogen, with carbon dioxide and ethane, and with carbon dioxide, nitrogen 
and ethane. 
3) To evaluate the available methods of prediction and to extend and 
improve the methods of Prediction. for Ene enthalpy of fluid mixtures. 
An isothermal throttling calorimeter was chosen for measuring 
the effect of pressure on enthalpy. In common with several other throttling 
calorimeters, a capillary tube is used to cause the pressure drop. An 
insulated Nichrome wire which passes inside the capillary for its entire 
length serves to supply electrical energy. This arrangement makes the 
expansion essentially isothermal. The pressure drop across the calorimeter 
is measured by a differential pressure transducer and any temperature 
difference between the inlet and outlet is indicated by a calibrated 
thermopile. 
Measurements were made with nitrogen as a test gas at four 
different temperatures from 273.15 to 374.15 K, and at pressures from 14 
to 140 bar. The data are in good agreement with values from the litera- 


ture. Measurements were made on the following six systems at pressures 
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from 14 to 140 bar: a) 14.5 mole percent of methane in carbon dioxide at 
temperatures of 0, 10, 20, 40, 60 and 90°C, b) 42.3 mole percent of 
methane in carbon dioxide at temperatures of 0, 40, 60 and 90°C, c) 
Equimolal mixture of methane, nitrogen and carbon dioxide at temperatures 
Onn On ancde Oa Caad ) Equimolal mixture of methane, carbon dioxide 

and ethane at temperatures of -10, 0, 20, 40, 60 and 90°C, e) methane rich 
mixture with carbon dioxide and ethane at Lenperatunes of s—10)5 0,520. 940. 
60 and 90°C and, f) Equimolal mixture of methane, carbon dioxide, nitrogen, 
and ethane at temperatures of -30, -10, 10. 30 and 60" GC, 

Two enthalpy tables and diagrams for the two mixtures of methane 
and carbon dioxide were constructed which are based on this experimental 
work and that of Bishnoi. 

The availability of more accurate enthalpy data based on 
calorimetric determinations provides an improved basis for comparison of 
available methods of prediction and the testing of the mixing rules. 

Nine enthalpy correlations were evaluated and compared. Three of the 
correlations used equations of state and six of the correlations used 

the corresponding states method of Curl and Pitzer with six different 
mixing rules. The six different mixing rules are: 1) Method of Barner 
and Quinlan, 2) Proposed Method, 3) Modified Method of Gunn, 4) Method 
of Prausnitz and Gunn, 5) Method of Joffe-Stewart, Burkhardt and Voo, and 
6) Method of Leland and Mueller. The three equations of state are 
Starling-BWR, Soave-Redlich-Kwong and Mark V correlation. It was found 
that the equation of state of Starling-BWR, Soave-Redlich-Kwong, Mark V 
and the corresponding states correlation with pseudocritical parameters 
estimated by method 2, are the most accurate of the nine enthalpy 


correlations tested. 
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INTRODUCTION 

A knowledge of the enthalpies of fluid mixtures over a wide 
range of pressure and temperature is necessary for engineering designs 
and calculations of thermal processes. In the past, enthalpies used in 
process design were obtained principally from P-V-T data, either directly 
using thermodynamic relations or indirectly using equations of state. 

The process of generating enthalpies from compressibility data involves 
differentiation with attendant loss in accuracy of at least one order of 
magnitude. In a region where the derivatives are changing rapidly, such 
as the critical region, the error may be very large. In addition, 
accurate volumetric data for mixtures is still rather scarce, and the 
determination of enthalpy changes across the two-phase region involves 
the use of not only volumetric data and their derivatives, but also 
vapor-liquid equilibrium data and derivatives. For this reason it is 
desirable to have direct experimental determinations of the enthalpy 
behavior of fluid mixtures under pressure. These data are quite scarce 
in the literature. 

The goal of this work was to construct a flow type calorimeter 
in a recycle system for the measurement of the isothermal effect of 
pressure on enthalpy. The isobaric effect of temperature on two 
binary mixtures of carbon dioxide and methane had been investigated by 
Bishnoi 4 . A specific purpose of the research has been to measure the 
isothermal effect of pressure on enthalpy of these two mixtures, so that 
it provides a thermodynamic consistency check and allows the construction 
of enthalpy tables for these two mixtures based on experimental data. In 
addition, experimental data of the isothermal effect of pressure on 


enthalpy were obtained for the following four systems: a) Equimolal carbon 
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dioxide-methane-nitrogen, (b) Equimolal carbon dioxide-methane-ethane, 
(c) methane-rich carbon dioxide-methane-ethane and (d) Equimolal carbon 
dioxide-methane-nitrogen-ethane. 

There are many different methods of prediction which have been 
proposed for enthalpies of fluid mixtures. A comparative study has been 
made by Natural Gas Processors Association’ ” on these numerous enthalpy 
correlations. The availability of more accurate enthalpy data based on 
calorimetric determinations provides an improved basis for comparison of 
available methods of prediction and the testing of the mixing rules. 

To reiterate, the specific goals of the present research were 
(1) to construct a calorimeter for obtaining the isothermal effect of 
pressure on enthalpy, (2) to make experimental enthalpy determinations 
on two carbon dioxide-methane mixtures, one carbon dioxide-methane- 
nitrogen mixture, two carbon dioxide-methane-ethane mixtures and one 
carbon dioxide-methane-nitrogen-ethane mixture, (3) to evaluate the 
available methods of prediction and (4) to extend and improve the method 


of prediction for the enthalpy of fluid mixtures. 
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OHCTION I - PRELIMINARY CONSIDERATIONS 
In this section, a discussion of flow calorimetry and the 
necessary equations for ¢ which are applied, and a discussion of » data 


in the literature are presented. 


Flow Calorimeters 

Measurements of the enthalpy of fluids at elevated pressures 
have been made by a number of methods. A review of experimental methods 
was presented by Barieau : - Flow calorimeters have many advantages and 
have been used widely. The calorimeter may be designed to operate in a 
number of differing modes depending upon the type of enthalpy data 
desired. 


The first law of thermodynamics, applied to a flow calorimeter 


with negligible potential and kinetic energy effects, is 


[Hp op ~ Hp ip J, 5 = (1) 
RIN) Vit 


where Q is the rate of the heat transfer, W the rate of electrical energy 
transfer, and F is the mass flow rate. In the isobaric mode of operation, 
the pressure difference P,-P, is made small and the fluid is heated to 


ail 


change its temperature. Equation (1) becomes: 


WI 2 
me ==- —- — —_— P 2 
En Hy te. x F | Gp | (2) 


where the integral term is a small correction for the fact that the 
pressure is not constant, and Q is assumed to be negligible. 


In the isenthalpic mode, no energy is added to the system and the 


heat leak is made negligible. For this case Equation (1) reduces to 


[H - H je ="0 (3) 
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In the isothermal mode of operation, a pressure drop is 
imposed on the fluid, and energy is added so that the outlet temperature 
is equal to the inlet temperature. It is practical to utilize this scheme 
only when the Joule-Thomson coefficient is positive, i.e. when the fluid 
cools upon expansion. For an isothermal operation, Equation (1) becomes: 
: ik 
eo i % Jt, x P be a ey 
i 
where 0 is assumed negligible and the integral corrects for any mismatch 
between the inlet and outlet temperatures. In this operation, a knowledge 
of the flow rate, the energy input and pressure drop are needed to determine 
the isothermal effect of pressure on enthalpy. 
The enthalpy change on mixing can be determined in a flow 
calorimeter by mixing two pure gases in a chamber and adding electrical 


energy to equalize the inlet and outlet temperatures. The first law of 


thermodynamics for this calorimeter reduces to: 


[H 


ae | [x,H, + x,H,] 
—mi1x 


Ty »P5 cell 22 T)>P, 


yg Sek (5) 
F 
Corrections can be made for the differences in pressure and temperature 


between the inlet and outlet. The excess enthalpy or heat of mixing can 


be determined by: 
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Thermodynante Relations 


The relations presented above involve integral changes in 
enthalpy. These integral data in a single phase region may be interpreted 


to yield the derivative properties. Point values of the isothermal 
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throttling coefficient $ are estimated from the integral values by 


application of the following relation from Equation (4), 


H — 
oa ce iy bs lim [ee] (7) 
iO vas eke apne 


The isothermal throttling coefficient can be estimated from 
known values of other thermodynamic properties. Thus, in making use of 


P-V-T data the following identity applies: 
b= v- 76, (8) 


The identity involving u, the Joule-Thomson coefficient and Ch» 
the isobaric heat capacity, provides another method of calculating 9. 


The functional relations between the enthalpy, the pressure and the temper- 


ature are: 


i ae 

ce b). aly Gals 5 it 
or 

0 ea ee 

GP) > y “oT p 
or 

a? 9 

o uC, (9) 

The virial equation of state is a power series in density: 

RT ava 2 sc ) 
The terms B, C, ... are known as the second, third, ... virial coefficients 


and they are related by the statistical theory of imperfect gases to the 


interactions of molecules in pairs, triplets, etc. 
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A similar power series is an expansion in pressure: 


Pye 2 
ee ee Decne ue pe (11) 


The coefficients of the two series are related by 


eae 
Stara 
2 
ct = €-~B) (12) 
(RT) 


Using Equation (8), (10), (11) and (12), 4 can be expressed as: 


2. dB! dc' 
¢ = -RT Cam rely pyar 3 She 
or 
dB ite ike Rn dc 
EE Oa arg sO Oe rag Cla Terral ey 


The zero-pressure value of the isothermal throttling coefficient 


is finite and depends only on B: 
¢° = B - TC; (14) 


Experimental values of $° can be used to derive changes in B upon 


integration of Equation (14): 


Ty F 
Be = T, | o°(t)dt + P| (15) 
Ty 


where t = 1/T. Experimental values of $° can be compared with values of 
6° calculated from virial coefficients obtained from other type of 


experiments. For mixtures, B is of the form: 


x oCeOe (16) 
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where Bay and ie are the pure component second virial coefficients and 


Be is the interaction second virial coefficient. It follows that the 


derivative of B. is: 
mix 
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Previous Expertmental Data for $ 


All of the above methods have been used to obtain thermal data 
of fluids at elevated pressures. Most of the results have been for pure 
components. A recent review of available pure components and mixtures 
data at elevated pressures have been presented by Mather 2 and Yesavage Hs 
Table I presents a listing of the experimental determinations of thermal 
properties under pressure by isothermal calorimeter, supplementing that 
of Mather a . Although an increased amount of enthalpy data are available, 


the need for more data to test methods of prediction and improve 


correlations remains. 
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TABLE I 


Recent References to Isothermal Effect of 


Pressure on Enthalpy of Fluid Under Pressure 


Vea System Author References 
1967 Nitrogen 

Propane-methane Mather et val. Sia iy OS 
1968 Propane 

Propane-methane Yesavage ef al. 9/,.98, 99. 100 
963 Methane-propane Dillard et ai. oe 
1969, 70 Carbon dioxide Vukalovich et al. 86, 87, 88 
L970 Methane-ethane Alkasab 2 
1970 Me thane-ethane-propane Furtado etal 26 
Loy Carbon dioxide-nitrogen Vukalovich et al. 89 
Leos Ethane Miyazaki 61a 
1974 Methane-carbon dioxide 


Methane-carbon dioxide-ethane 


Carbon dioxide Peterson & Wilson 68 
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SECTION II - THE CALORIMETER AND THE FLOW SYSTEM 
A review of the previous work employing isothermal throttling 
calorimeters is presented here. The design of a new throttling calorimeter 


and the experimental equipment are discussed in this section. 


Baekground on Isothermal Throttling Calorimetry 


All methods for measuring the isothermal effect of pressure on 
enthalpy involve some device to cause the pressure drop and a source of 
energy to compensate for the cooling effect resulting from the pressure 
drop. So, the methods are suitable only when o is negative and no 
isothermal calorimeters have been designed to extract energy if the 
fluid warms on expansion. 

A review of the development in isothermal throttling calorimetry 
had been presented by Mather =a - In general, a porous media, a needle 
valve or a capillary tube is used to cause the pressure drop, and the 
electric heating coil is used as a source of energy. 

The calorimeter designed by Mather Zt consists of a capillary 
coil to cause the pressure drop and the insulated heating wire is placed 
inside the capillary. Measurements were made on nitrogen and a mixture 
of propane and methane in the temperature range from -140°F to 200°F and 
in the pressure range from 100 to 2000 psia. Yesavage i used the same 
equipment to make measurements on propane and propane-methane mixtures. 

Vukalovich Ea aes reported the measurements on carbon 
dioxide by using adjustable needle valve with pressure drop up to 50 atm. 
Later on, TUccat yicheetpaleon designed a different type of calorimeter 
in which the throttling unit is pressure-relieved. Again the adjustable 


needle valve was used as the throttling unit. 
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Alkasab designed an adjustable needle valve to study the 
isothermal throttling effect of methane-ethane mixtures. 

: 68 ; : ; 

Peterson and Wilson constructed an integral isothermal throttling 
calorimeter by using throttling valve to cause the pressure drop. Measure- 
ments were made on pure carbon dioxide, and mixtures of equimolal methane- 
carbon dioxide, equimolal methane-ethane-carbon dioxide. 

Fj 25 ‘ ; ‘ ; 

Eakin and DeVaney designed a calorimeter which permits 
measurement of the Joule-Thomson coefficient, the isothermal enthalpy 
change on throttling and the isobaric enthalpy change. The idea of 
tubing packed with powder was used to cause the pressure drop. This 


idea was first used for Joule-Thomson experiments by Dawe and sue omae 


Destgn of Isothermal Throttling Calortmeter 


The first requirement in the design of the calorimeter was that 
it should operate over a wide range of temperatures from -100°C to 200°C 
and at any pressure in this temperature range fromel00 tor 2000 %psia. 
Stainless steel was selected as the material of construction because of 
anticipated subsequent experimental work with corrosive gases and its 
strength at low temperature. 

Since the design was to allow for the possibility of measuring 
isothermal enthalpy changes across a two-phase region, the use of throt- 
tling valve was abandoned because the instability of two-phase flow 
through a valve caused by alternate slugs of liquid and vapor prevents 
accurate thermodynamic measurements. 

A capillary was considered to be the best configuration for the 
throttle because the expansion is smooth and if energy is added to the 


fluid during expansion, the process will be almost isothermal. Since 
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ET 
enthalpy is a point function, it is not necessary for the expansion to be 
isothermal. However, by making the path as nearly isothermal as possible, 
heat leakage from the Surroundings to the system is minimized. 

The disadvantage of capillary throttle is that it is not 
possible to vary flow rate and pressure drop independently. The disadvan- 
tage can be overcome by making the capillary removable so that various 
combinations of pressure drop and flow rate may be obtained. 

The heater of the present calorimeter is an insulated Nichrome 
wire which passes inside the capillary for its entire length. Two copper 
power leads are connected to the heater wire at the ends of the capillary. 
In order to minimize the power losses through the leads the current must 
be kept low. The combination of heater size and capillary size allows 
different pressure drop and power input relations. With this arrangement, 


currents were kept below two amperes. 


Desertptton of the Calorimeter 


The isothermal throttling calorimeter of this work is shown in 
Figure 1. The fluid enters the vacuum jacket from the constant temperature 
bath. The pressure tap at the inlet portion of the calorimeter serves to 
measure the inlet pressure and a differential thermocouple is used to check 
the temperature between the inlet portion of the calorimeter and the bath. 
Then the fluid passes through the throttling portion. A removable 
capillary coil at the throttling portion causes the pressure drop. The 
outlet temperature of the fluid is brought back to the inlet temperature 
by applying power to the Nichrome wire, which is placed inside the 
capillary. In this work, three different capillary sizes are employed of 


15 BWG, 16 BWG and 17 BWG hypodermic tubing, all about ten feet in 
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Assembly Drawing for Isothermal Throttling Calorimeter 


FIGURE 1. 
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113 
length. Two different sizes of Nichrome wire, 36 and 40 gauge, are used. 
A five junction differential thermopile is used to determine the equality 
of inlet and outlet temperature. The pressure drop is determined by the 
differential pressure transducer. As shown in Figure 1, all the 
pressure and temperature measurements are at the same elevation. The 
lead wires of the heater and the thermopile are wrapped around a massive 
aluminum thermal equalization block. The heater wire leads are brought 
from the vacuum jacket into both ends of the capillary by Conax connectors. 
The detail drawing of this portion is shown in Figure 2. 

Other efforts were made to reduce heat transfer to the surround- 
ings. A radiation shield completely enclosed the Capillary coil. The 
power leads and thermopile leads are brought into contact with the 
vacuum jacket. The entire jacket is evacuated to less than five microns 


through the vacuum line. 


Eapertmental Equipment 


The calorimeter is part of a recycle system, the evolution of 
which can be traced in the thesis of Bishnoi 2 . The recycle system 
serves to bring the fluid under investigation to the desired conditions 
of pressure and temperature for measurements, A schematic diagram of the 
flow system is shown in Figure 3 . The fluid delivered by the compressor 
goes to the high pressure storage tanks, one of which is connected to the 
low pressure surge tank through the by-pass control valve BPCY. The 
fluid leaving the surge tanks goes to a hand valve and the high pressure 
control valve HPCV. The fluid passes through a long cooling coil before 
it enters the calorimeter. The cooling coil and the calorimeter are in 


the same bath. The pressure at the inlet of the calorimeter is governed 
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. Inlet tubing for the Calorimeter 
. Isothermal Calorimeter 
. Capillary Tubing 

. Liquid Nitrogen Cooling 
. Controlled Heat Input 


Stirrer 


. Platinum Thermometer 
. Sensor for Heat Input Controller 

. Temperature Bath for the Calorimeter 
. Joule Thomson Control Valve 
. Water Cooling 


. Hand Valve 


. Orifice Meter 

. Temperature Bath for the Flow Meter 
. Low Pressure Tank 

. Storage Tank 

. Corblin Diaphragm Compressor 

. By-pass Control Valve 

. High Pressure Tank 

. High Pressure Control Valve 

. Feed Valve 

. Discharge Valve 


FIGURE 3. Flow Diagram of the Apparatus 
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by the set point of the high pressure controller. The outlet stream of 

the calorimeter goes through a hand valve and the control valve JTCV. The 
opening of JTCV, manipulated manually, gives the desired mass flow rate 

of the fluid through the calorimeter. The fluid is then passed through 

the flow meter. The constant inlet pressure of the flow meter is controlled 
by the BPCV valve which is governed by the set point of the controller. 

Then the fluid passes through a low Pressure storage tank before returning 


to the intake of the compressor. 


Auxtltary Equipment 


Z. Control Valves and Pressure Control 

It was necessary to control the pressure of fluid entering the 
calorimeter and the flow meter. For flexibility, the low pressure control 
scheme should be such that it could maintain the pressure anywhere in the 
range of 10 to 50 psig. The high pressure control scheme was required 
to maintain any pressure up to about 3000 psig. Both the schemes should 
also be able to handle flow rates ranging from one to four standard cubic 
feet per minute for different gases such as methane, nitrogen, carbon 
dioxide, ethane or their mixtures. 

Figure 3 gives the location of HPCV, the high pressure control 
valve, JTCV, the pressure reducing control valve and BPCV, the by-pass 
control valve. These valves were sized by the standard method of Sy 


(valve flow coefficient) calculations. The valves used are: 


HPCV - Annin valve model 9460, trim size 0.02, air-to-close, can 
handle fluids from -400 to 750°F. 
eG - Annin valve model 9460, trim size 0.003, air-to-open, 


can handle fluids from =400 to 4/50°F. 
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BPCV -  Annin valve model 5061 "Wee-Willie", trim size "A" - C. 
range 0.0001 to 0.01, field reversible, can handle fluids 
frome=200 > tol465 °F. 

The valves had pneumatic positioners and stainless steel type 
316 bodies. 

For both HPCV and BPCYV, single feed back control loops having a 
controller with a proportional constant and an integral constant were 
used. Foxboro controllers with built-in power supplies were used for the 
high and low pressure controls. High and low pressures were measured, at 
the locations shown in Figure 3 , by Foxboro pressure transmitters model 
611 GH and 611 GM respectively. The current to pressure converters used 
in the pressure control loop are Fisher type 546 electro-pneumatic trans- 
ducers. 

The resistance of the valve JTCV was adjusted to give the 
desired flow rate by manually setting a potentiometer which supplied a 
current signal to the current-to-pressure conyerter, giving an air 
pressure to the valve positioner. Since the upstream pressure of JTCV 
and its resistance remained constant and the flow through the valve was 


sonic, this scheme gives constant mass flow rate through the calorimeter. 


2. Matertal Rectrculating Device 

A two-stage diaphragm compressor was used to develop the 
necessary pressure and to circulate the fluids. It was a Corblin A2CV250 
which was rated for a discharge pressure of up to 3700 psig with a flow 
capacity of 3.5 scfm at a suction pressure of 14.7 psia. The inlet and 
delivery lines of the compressor were connected to the tanks by means of 


flexible pipes. 
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5. Gas Heaters 

In order to maintain the steady-state operation, the control 
valves HPCV, BPCV and JTCV were maintained at constant opening. The 
pressure drop across the "HPCV" valve ranged from 50 to 400 psi and a strip 
heater was used to compensate the cooling effect due to the pressure 
drop. Since cooling effect is large because of pressure drop across the 
BPCV and JTCV ranged from 300 to 2000 psi, steam heaters were used to 


heat the fluid before it passed through the valves. 


4. Temperature Controller 
The calorimeter bath temperatures were controlled within +0.01°C 


by a Hallikainen Thermotrol temperature controller, model 1053. 


Measuring Instruments 


1. The temperature of the bath was determined by using a Leeds and 
Northrup platinum resistance thermometer, model number 8163 and it was 
calibrated for use in the range 90.188 K < T < 773.15 K on the International 
practical temperature scale of 1968 by National Research Council of 
Canada. 

2. The inlet temperature of the fluid to the calorimeter was assumed to 
be equal to the bath temperature and a one-junction copper-constantan 
differential thermocouple was used to check the measurement. The 
temperature difference between inlet and outlet was measured by a 
five-junction copper-constantan differential thermopile. Both the 

one- junction thermocouple and the five-junction thermopile were calibrated 
for the zero reading at boiling and ice temperatures of water, and the 


melting temperature of dry ice. A d-c null voltmeter was used to read 
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19 
the null point for the thermopile. The accuracy of the temperature 
measurement is about +0.01°C. 

3. The inlet pressure was measured by a Heise gauge and the pressure drop 
across the calorimeter was determined by a Validyne differential pressure 
transducer. Both Heise gauge and pressure transducer were calibrated 

by the Ruska dead weight tester. The accuracy of the gauge is 0.1 
percent, of full scale, and the accuracy of the pressure transducer is 0.5 
percent of the full, scale, 

4. The electrical energy input to the calorimeter was supplied by Kepco 
DC power supply and was measured by a Hewlett-Packard model 3450 
multifunction meter with the six-digit display. 

5. In the early stage of this work a Nupro "F" inline filter was used 

as a flow meter. The mass flow rate for nitrogen and 42.3% methane in 
carbon dioxide mixture was determined from the inlet pressure and the 
pressure drop across the filter together with the constant temperature at 


the element. These data were used to solve the calibration equation: 


2 2 
Peo b 
dl 2 FM 
soe Sa i 8 


for the mass flow rate F. The constants A and B were obtained from a 
least square fit of the calibration data obtained by the gasometer, 
where P, and P, are the inlet and outlet pressure of the flow meter which 
are given by pressure transducers. wu', Z and M are the viscosity, 
compressibility and molecular weight of the gas respectively. 

During the data taking of the 42.3% methane in carbon dioxide 
mixture, it was found that the solid contamination in the flow meter 


element causes the changing of the calibration curve. The flow meter 


was recalibrated and after obtaining the experimental data of the mixture, 
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20 
the flow meter was recalibrated again. It was found that the two calibra- 
tion curves disagreed by 0.8 percent. So, the filter was abandoned as 
device for the flow measurement. 

A square-edged orifice was then used for the flow measurement. 
The orifice had a diameter of 0.16 in. and it was located in a laine 
tube with an inside diameter of 0.429 in. The inlet pressure was 
maintained constant by the BPCV control valve, and was determined by a 
Validyne pressure transducer. The pressure drop across the orifice was 
measured by a Foxboro d/p cell transmitter type 613 M. The accuracy of 
the pressure transducer and d/p cell is 0.5 percent of the full scale. 
The orifice was calibrated using pure nitrogen, methane, ethane and carbon 
dioxide and argon by the gasometer. The orifice was calibrated at the 
completion of the experiments on every mixture using the gas mixture 
itself. The overall standard deviation of the flow rate is about 0.5 
percent. The calibration of the orifice is given in Appendix A. 

All the output signals of the pressure transducers, d/p cell 
transmitter, platinum resistance thermometer, power to the calorimeter 
and thermopile were connected to a Hewlett-Packard model 3450 multi- 
function meter by a multi-junction selection switch. 

Figure 4 shows the electric circuits of platinum resistance 


thermometer and the calorimeter heater. 
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Figure 4: Wiring Diagram of Platinum Resistance Thermometer and 


Calorimeter Heater Circuits 
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SECTION III - EXPERIMENTAL RESULTS AND DISCUSSION 
In this section, the experimental procedure, the materials used, 
and a summary of the techniques used in interpreting the experimental 
data obtained in this investigation are presented. The experimental 
results of » of nitrogen, two mixtures of methane and carbon dioxide, two 
mixtures of methane, carbon dioxide and ethane, one mixture of methane, 
carbon dioxide and nitrogen, and a mixture of methane, carbon dioxide, 


nitrogen and ethane at various temperatures are also presented. 


Procedure of the Experimental Measurements 


The fluid under study was brought to the desired state of 
pressure and temperature at the inlet to the calorimeter with the recycle 
system shown in Figure 3. A pressure drop of 60 to 250 psi was obtained 
by adjustment of the flow of fluid through the calorimeter section. The 
flow rate was adjusted by control yalve JICV, manipulated manually. The 
voltage supplied to the calorimeter heater wire was adjusted manually 
until the temperature difference between the inlet and outlet of the 
calorimeter, indicated by the differential thermopiles was less than two 
microvolts (0.01°C). When steady state was reached, as indicated by 
the constancy of the pressure drop and the zero temperature difference 
across the calorimeter, the values of the variables were recorded. The 


time of approach to steady state was of the order of half to one hour. 


Materials Used 


The source and purity of the gases used in this work are given 


in Table II. The gases were used without further purification. 
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Component 


Nitrogen 
Methane 
Carbon dioxide 


Ethane 


TABLE. IT 


Materials 


Supplier 


Consumer's Welding Co. 


Matheson Gas Products 


Canadian Liquid Air Co. 


Matheson Gas Products 


Ltd. 
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Purity 


99.9932 
OD elaine 
99.95 % 


Oo 90 
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Composttton Determitnatton 


The composition of the mixtures in this study was determined 
chromatographically. The calibrations for the gas chromatograph are 
presented in Appendix A. The precision of the composition determinations 
was approximately +0.2 mole percent as was checked by making several runs 


for each isotherm. 


Interpretatton of Data 


The basic data were recorded in terms of quantities that can 
be readily measured, such as microvolts and height of a fluid. These 
quantities were converted to temperature, pressure, pressure drop, power 
input and flow rate. In the single phase region, the isothermal throttling 
coefficient » which is shown in Equation (7) is obtained from the 
integral data by several techniques. Both graphical and computer reduction 


methods were used. 


Grapht cal Interpretation 


One may consider determining » at constant temperature as a 
function of pressure. The experimental data were reported as sets of T, 
P P, an cs over a wide range of temperature. For any one 


data point, 
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where or, is the mean isothermal throttling coefficient between Pi and P,: 


Values of On were plotted as horizontal lines are shown in Figures 5 and 
m 


10. The smooth bn curve was obtained by satisfying Equation (19): the 
m 


area under the horizontal line segments, on *(P, 
m 


area under the smooth isothermal throttling coefficient curve. 


= Pi)> should equal the 


Computer Reduction 


The above graphical procedures are extremely time consuming 
and in addition the equal area construction can easily lead to errors. 
Therefore, a computer program was developed for interpretation of integral 


data. The true isothermal throttling coefficient is given by, 


lim qa oH 
aa latina Ube ian EET 


where q = His the rate of electrical energy transfer per unit mass flow 
rate. The result of actual measurements is the mean Oe Om ae, = Pi)» 
associated with the mean pressure of the interval, haba (PS ae P,)/2. 
Deviation from linearity of the ¢ versus P curve will therefore require 
adjustment of the mean isothermal throttling coefficient by a curvature 
correction to yield the true isothermal throttling coefficient at ideae 

The curvature correction can be determined by considering the 
isothermal throttling coefficient to be represented by a cubic equation 


in pressure over a limited range 
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OEM. pees hae 3 ee) oe 
= a) (P, - PD Vee oats P,) + Aes eae ig eS ee eee) 
and the mean ¢ is, 
*m PP 
2 1 
a a 
= Santee oy 2 
ay aE 5) (P, 5 0 P,) ae 3 (Py aR PoPL =F PD 
eect ew) Dies okt 
es Cicy byl ee (23) 
The difference between the true ¢ and the mean dm at Py is, 
see pe 
Oras Olan cr [2a, ats 3a,(P, SL Pi] (24) 
Since at the mean pressure, 
o 
Ces ae ee (ere (25) 
2 Z Siam: 2 
oP 
the true isothermal throttling coefficient is, 
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ee) aman 
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values of ( 


yr were obtained by considering the mean isothermal to be 


oP 
sufficiently close to the true isothermal throttling coefficient and 
differentiating the smooth mean isothermal throttling coefficient. 
The extended spline fit technique was used to obtain a smooth 


mean isothermal throttling coefficient and its second derivative with 


Lespect to Pp. 


The Extended Spline Fit Method 


The spline fit technique has been discussed by Landis and 


Nilson as - This method puts a different cubic between every two 
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2 
successive data points such that the curve passes exactly through each 
data point and that the first two derivatives of the curve on the right 
hand side of the data point are equal, respectively, to the first two 
derivatives of the curve on the left hand Side of the data point, all 
derivatives evaluated at the data point. 

If, instead of defining the interval boundaries to pass 
through every data point, interval boundaries are determined arbitrarily 
such that each interval may contain a number of data points, the method 
can be extended to provide smoothing of the data. The details of this 


method are discussed by Klaus and Van Ness ae : 


katenston to Low Pressure 

Since the lower limit in pressure of the recycle flow system 
is 200 psia, data from the literature are used in extending the enthalpy 
results to zero pressure. The experimental isothermal throttling 
coefficients obtained at elevated pressures should extrapolate to the 
zero pressure values derived from the experimental second virial 
coefficients using Equation (14). The resulting curve was integrated to 


determine the effect of pressure on enthalpy at low pressures. 
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Expertmental Results With Nitrogen 


The experimental data were taken for nitrogen on the isothermal 
throttling calorimeter at four different temperatures from 273.15 to 
374.15 K, and the pressure up to 2400 psia. Figure 5 shows the isothermal 
throttling coefficient, » at 273.15 K plotted as a function of pressure. 
An equal area curve passing through the bars obtained point values of ¢. 
The agreement between this work and Maeheroe is less than 1.0%. Figure 
6 and 7 show the comparison of two different temperatures with other 
experimental results. The experimental band indicates the precision of 
about +1/2% for this experiment. The data of Ishkin and herman” are 
believed to be high by as much as 5% (see Verdes) 

The values of $ can be calculated from an equation of state 
using Equation (8). The Benedict-Webb-Rubin® equation is frequently used 
for calculation of thermodynamic properties. For nitrogen, two sets of 
constants for this equation are available, those of Stotler and Renedwer 
and those of Bloomer and mag: Isothermal throttling coefficients were 
also calculated from the equation of state for nitrogen from the U.S. 
Bureau of ear The results of equation of state calculation were 
compared with the experimental data in Figure 8 and 9. 

Table III presents the smoothed experimental measurements of fo) 
by this work. 

Enthalpy departures were calculated from the experimental data 
by integration of » with respect to pressure. The results for the three 
isotherms of this work were compared with other experimental calorimetric 


results in Table IV. The present data were in very good agreement with 


the most recent experimental values in the literature. 
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TABLE III 


ISOTHERMAL THROTTLING COEFFICIENTS FOR 
NITROGEN 


=p sail © bee 


Temperature, °C 


PRESSURE 

(bar) 0 30 70.3 HORE 
O* 7.48 Bloe)s) 4.58 5,02 

10 ees 5.83 4.44 3.51 
20 Teale. Oo 4.30 3.40 
30 1203 5.54 4.16 eas) 
40 6.86 Diedd 4.01 Shiny hy) 
50 6.66 evoyedl 387 3.04 
60 6.46 5.04 Se! Pies 
70 6.26 4.89 3.98 Zeod 
80 6.06 4.73 3.43 Post [i 
90 5.85 4.55 Seal 2.60 
100 5.63 4.37 3.18 2.48 
110 559 4.20 3.05 EC CoW | 
120 5.14 4.01 Pa ESPNS 
130 4.88 3.83 Pee 2.14 
140 4.75 320) 2.6/7 2.04 
LoO 4.34 3.46 4 ey) To 
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TABLE LV 
COMPARISON OF EXPERIMENTAL ENTHALPY DEPARTURES 


FOR NITROGEN WITH OTHER EXPERIMENTAL WORKS 


MCS WS Ce) Ghee ey Sie TSU LES in Cie Gea eso 
(H.-H TF : J mol : 
TEMPERATURE PRESSURE P 


(K) (bar) This Work ethers Dawe & Suewdenas 

Pad) Byles) 20 146.6 147 148 
40 28740 290 291 

60 420.2 422 428 

80 545.4 543 558 

100 662.4 660 679 

si0s hale) 20 L1LGwa - ny) 
40 2272 - 230 

60 set eRS) = 338 

80 429 32 = 439 

100 o2 ez = oyeV 

3/4. 15 20 V0.2 - 67 
40 136.0 = 132 

60 196.8 - 194 

80 25 GeL = PASH! 


100 SEG) = 305 


36 


hxpertmental Results With Methane - Carbon Dtioxtde Mixtures 


The measurements of the isothermal throttling coefficient were 
made on 42.3 mole percent of methane in carbon dioxide at four isotherms, 
Oye O en o0n and. 9020.6. The pressure ranged from 400 to 2000 psia. Measure- 
ments were also made at six isotherms and the same pressure range for 
14.5 mole percent of methane in carbon dioxide. The six isotherms were 
OS 10, 205 #40. GO and 90°C. 

Typical isothermal data are presented in Figure 10. Average 
values of om = (AH/AP),, were plotted as horizontal lines and the computer 
reduction method was used to determine point values of » = f(P) as 
illustrated by the solid curve. In extrapolating the results of the 
experimental investigation to zero pressure, it was necessary to use 
Equation (14) for $° in terms of second virial coefficient. Values of o° 
were calculated from the virial coefficients given by Dymond and Snare 
and Neo The calculation of $¢° is shown in Appendix C. 

Tables of ¢ values of all of the experimental isotherms are 


presented in Table V and VI. Figure 11 and 12 present the plot of ¢ as a 


function of temperature at different isobars for the two mixtures. 


Compartson of Results 


Since the isothermal experimental data on the carbon dioxide- 
methane mixtures with same compositions are not available in the literature, 
comparison cannot be made. However, the isothermal throttling coefficients 
of two mixtures at each isotherm and at a pressure up to 2000 psia can be 
calculated by using BWR equation of state and Bishnoi and Robinson? 
mixing rules. The parameters of carbon dioxide and methane, reported by 


Bishnoi and Ropineon.— were utilized in the above calculations. 
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FIGURE 10. Isothermal Throttling Coefficient for the 14.5 percent CH, 
in CO5 Mixture at 40°C 
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TABLE V 
ISOTHERMAL THROTTLING COEFFICIENTS FOR 
THE 42.3 MOLE PERCENT METHANE IN CARBON DIOXIDE 


-¢ (J Asie fate 


Temperature °C 


PRESSURE 0) 40 60 90 
(bar) 

Ox 32246 24.10 Zea Teo 
10 34.42 24.95 21.92 18.48 
70) S/o O 269) 22.58 18.78 
30 40.70 Die. 09 23520 19),,05 
40 45.41 28 = 27 23.95 Tora 
50 BA iss 29.48 24.60 19.58 
60 62.04 30.76 25.20 19.83 
70 al 32.05 25.76 207,01 
80 8/615 33.24 Oe 20).13 
90 biog! 34.23 26.69 20. 107, 

100 iad! 34.90 26.95 20517, 
110 43.52 34.97 27.00 20 42 
120 30.04 34.00 26.76 20.00 
130 Zio 5 32.27 2601 /, 19.74 
137 ore oleng Zo eo3 ia3 
"po pee 


TABLE VI 
ISOTHERMAL THROTTLING COEFFICIENTS FOR 
THE 14.5 MOLE PERCENT METHANE IN CARBON DIOXIDE 


=¢ (J not bare} 


Temperature, °C 


PRESSURE 0 10 20 40 60 
(bar) 

O* 44.23 40.92 50679 eyes 27.78 
10 49.61 45.64 395.09 33.48 295 
20 Sid) MS 51-380 43.82 5) 30.67 
30 69.86 60.15 48.57 36.28 32.19 
40 = 1oege 54.82 42013 Sha Sh 
50 = OGn/3 67.63 46.32 36.04 
60 = = 90.03 51.54 38.37 
70 = = 12925 58.39 40.92 
80 og a 241.8 67.48 43.66 
90 TE SEM 20.90 99.40 78.06 46.53 

100 6.10 14.21 44.51 85.15 49.17 
Hai) 5.02 VOR 7S ie ifs 80.53 50.91 
120 3298 8.43 796 63.06 07, 
130 = 7.06 14.65 46.30 eel 
137 - 6.38 11.54 37.00 46.37 
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FIGURE 11. Isothermal Throttling Coefficient for the 42.3 percent CHy 
in COs Mixture at Various Constant Pressures 
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FIGURE 12. Isothermal Throttling Coefficient for the 14.5 percent CHg in CO5 
Mixture at Various Constant Pressures 
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In general, the deviations are larger at the lower pressures of 
the lower temperature. The comparison of the experimental data and the 
calculated values from BWR equation of state for 42.3 percent methane in 
carbon dioxide at 0°C and 14.5% mixture at 90°C are shown in Figure 15 
and 14 respectively. 

Peterson and ‘ilies reported results of isothermal enthalpy 
departure from isothermal throttling experiments for equimolal methane- 
carbon dioxide mixture at pressure from 20 to 2000 psia and temperature 
from -50°F to 300°F. Isothermal enthalpy departures of the two mixtures 
of this work were calculated by integration with respect to pressure. 

The results were then interpolated to 100°F and 150°F at pressure of 800, 
1200, 1600 and 2000 psia. The final results were plotted at constant 
pressure and temperature as a function of mole fraction. Enthalpy 
departure of pure methane and carbon dioxide were obtained from modified 


rath with the parameters from 


BWR equation proposed by Starling (BWRS) 
heen Figure 15 shows the comparison of the results between this 


work and that of Peterson and moleen | at 100°F. The comparison at 150°F 


is also in good agreement. 


Enthalpy Tables and Diagrams 


The experimental data on the effect of pressure on enthalpy by 
this work and the effect of temperature on enthalpy by Bienne a for the 
methane - carbon dioxide mixtures have been used to prepare skeleton 
enthalpy tables and diagrams. Tables and diagrams are presented for 14.5 
and 42.3 mole percent methane in carbon dioxide mixtures. Comparisons 


were made with excess enthalpy data from the literature. 


: | Ae —-_ 
sexinsiase, inductions RG) Sotto papi 
sunisteems) ban’ 6ieq oor: o7 0S mot smuessry a s nat 
. Bendeka awit Gy to astuiteqeb Ne, ir 3° 008 


jocguaot Te battol4 s9oW aripeny Tenet, Digs ening 08 a 
yqisaies .nottoos? eine Te neegonva ges Srusjsresanegt 
bolisham most banksido Ee abrxoth Adres hae so Sot 
sox? eyaienbisg Set silbe BV. 98 (oat) pe %d beaiqongit 
fai aesvosd ata tl eda Tm, nse lisa sia apie 2! sxpgit- _ 
B°0Cs 32 ruelreqego salt JAP a0 ry ar poke bax nopzsisT is ver bam ee 


em | ;Semaorye hee et - 


qu yqltains os ayoReezy to josiiy add uo sinh Jevasettsgxs sit 
ant 703 © ciliata vd. yglerians. co aaron to ‘tosis sat bas stow whds 
nossiste ‘stagexq 03 fem age oved eorwixte ahtxoih poives + samdaom 
-eostest’ bre asides yqhedaae. 
shadiem Ins eq alom €,Sé bas 


ab sii 2eeox9 tiv shes ssew 


43 


Ovi 


> GL'ELTZ Ye B4NI XII) 1U9d19q BOW EZ 40% “EL AYNDIS 


410q_ "FYNSSIYd 
ora 001 og 09 Ov 02 


ip 
ap 1-d = Vv 


ppiuesuiqoy pue touysig) YMg eS 


YOM SIYL 


LLG'0 apixoiq uoqieg uoNdes 4 BJO 
EZp'O suey UONDe14 BjOW 
GLELZ = 1 


i \noridoR bis onde) AWE ® 
96 yg = % A 
Th 


00! 08 os ob 
yod ,aRUecsR9 


St ENS 25 awixiM theswePeloMESdio® £f 3AUAIS 


os 


44 


> GL'E9E 3€ B4NIX IA) }UBII9Yg BIOW GyLjJO% “pL FHNOIS 


40q_*3Y9NSSIYd 
Ovi Oz1 001 08 09 7 


ip 
gp i-—4d = fo 6 


ypuosuidod pue fouysig) YMA @ 


HOM SIYL 


02 


GG8'0 ap!xoig uoqueg UONDe14 ajo 
GrL'O dueY IA] UOIJDCIY BjO|\j 


A GLE9E = 1 


pen cs eee wt 


x 


~~ 
ee _ 


- Ohr—ts«éST oor 08 08 On “o.o 
_ i, vod «=, 3ay2e555 
a M ZF. 806 ts euixiMtosmeteHM G40 toy OF BAU 


CHy — CO, 


3900 100°F 
® sBwrs2®? 


A Peterson & Wilson ©8 


3300 O This Work 

2000 psia 

2700 1600 psia 

1200 psia 

800 psia 
2100 
1500 
900 
300 


0 Of 0.4 0.6 0.8 1.0 
MOLE FRACTION OF CH, 


FIGURE 15. Enthalpy Departures for Mixtures of Methane 


and Carbon Dioxide 
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Nominal 14.5 Percent Methane in Carbon Dioxide Mixture 
A skeleton enthalpy table and a pressure-enthalpy-temperature 
diagram for this mixture have been prepared which are based entirely on 
experimental calorimetric data. The following procedure was used in 


preparing the diagram and table: 


b. Reference state of 0 J/mole were taken to be the pure components 


as perfect crystals under zero pressure at absolute zero temperature. 


Zn The enthalpy of pure methane as a gas at zero pressure and 0 K 
was calculated using the enthalpy of methane (saturated vapor) at 
100 K reported by igsie the B-W-R equation of state with parameters 
from Bishnoi and Repinaanes to correct from 0.344 atm to zero 

61 


pressure, and value of the ideal gas enthalpies from McBride etal. 


the following results were obtained: 


kJ moi? 
Enthalpy of methane at 100 K and 
0.344 atm i eeae To 
Corrected to ideal gas (0.344 to 0 
atm) on MEWS) 
Ideal gas enthalpy change from 100 
COnUeK Sa eX 
9.22025 


so, the enthalpy of pure methane from perfect crystal at 0 K to any 


temperature at zero pressure is, 


-1 
b— 9 _— © 
i, 9.22025 + (Hy, H) kJ mol 


where (Hy, = He) is the ideal gas enthalpy change reported by 


b 
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47 
McBride et a1°?, 


The enthalpy of pure carbon dioxide as a gas at zero pressure and 


any temperature is reported by Vukalovich and atermeltis as follows: 


a oly, oe ° =]L 
H = 26.2235 + (Hy Hq’) kJ mol 


where (H, - H) is the ideal gas enthalpy change reported by 
McBride et alo. 

The enthalpy of methane-carbon dioxide mixture at zero pressure 
and any temperature T K was calculated assuming zero heat of 
mixing under these conditions. 

Hy mix ~ : gered 


ke} cules 


23.7580 + 0.145 (A, = HS cH ots 0.855 (H., = Ho) 
4 


cO., 

The isothermal effect of pressure on enthalpy at) 293.15, 
313.15, 333.15 and 363.15 K was obtained from the integral > data 
of this work from 0 to 138 bar. For temperatures at 283.15 and 
273.15 K, the isothermal effect of pressure on enthalpy was obtained 
from 0 to 35 bar. The enthalpy departures for this mixture are 
presented as function of temperature in Figure 16 and comparison is 
made with values calculated from equations of state. BWR equation 
with Bishnoi and Ropineonu. mixing rules predicts higher departures, 
while Soave-Redlich-Kwone’ © equation with the interaction constant 
ko i. com aes is seen to predict lower departures. BWRS equation 
with parameters of Hopke and iy a gives better predictions. 


The isobaric heat capacity data of Stoner At ase OO e Ooo ke Ll, 


363.6 and 423.75 K was used to determine the isobaric effect of 
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48 
temperature on enthalpy. Since enthalpy is a property, changes in 
enthalpy are independent of the path chosen and the consistency of 
the experimental data can be tested by making loop checks. A grid 
consisting of 8 loops which is shown in Figure 34 was constructed 
and the enthalpy change around each loop evaluated. The actual 
sum divided by the sum of the absolute values of the enthalpy 


differences in the loop, 


% AHN 
a 


% deviation = =n. x 100 (27) 
Sect 
al 


provides a good measure of the consistency of the eae The 
maximum percentage deviation of each loop is less than 0.37 percent 


and the percentage deviation of each loop is shown in Table XXVII. 


Adjustments were made on individual values of AH. as 
required to make y AH. = 0 for all loops. These adjustments were 
made within the Ree ce of precision of the basic data. The 
accuracy of the heat capacity data is within AN yee Error 
analysis of this work is presented in Appendix A. The accuracy of 
the isothermal throttling coefficient is within +1.0%. The amount 
of the adjustments are presented in Figure 34. 


Ge A smooth plot of the results was prepared and is shown in 


Figure 17. The values of enthalpy are given in Table VII. 


Nominal 42.3 Percent Methane tn Carbon Dtoxtde Mixture 


A skeleton enthalpy table and a pressure-enthalpy diagram for 


this mixture have been prepared. 
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FIGURE 16. Enthalpy Departures for the 14.5 Percent Mixture 
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The isothermal effect of pressure on enGhalloygat 273215, 9313.15 
and 333.15K was obtained from the integrated » data of this work from 0 
jefe) ALSoy [ayebers 

The isobaric heat capacity data of Biennial | ate3l2 07,333.63 
and 363.45K was used to determine the isobaric effect of temperature on 
enthalpy. 

A grid consisting of 8 loops was constructed and enthalpy change 
around each loop evaluated. Figure 35 shows the location and enthalpy 
changes of the loops. In all closure the percentage deviation is less 
than 0.3 percent and the percentage deviation of each loop is shown in 
Table XXVIII. 

Adjustments were made to make : AH. ="O"fore ali oops™ and tne 
amount of adjustment is presented in oe Bs 

A smooth plot of the results was prepared and is shown in Figure 


18. The values of enthalpy are given in Table VIII. 


Comparison of Methane - Carbon DMoxtde Mixtures Enthalpy Data 


The heat of mixing or excess enthalpy for the two mixtures was 
calculated at various pressures from the data of Table VII and VIII 
together with the methane data of pene and carbon dioxide data of 
Vukolovich et al°>. The results were compared with the excess enthalpy 
data of Lee and MelcHerac directly determined using a flow calorimeter. 


The comparison is shown in Figure 19 and 20. The agreement is generally 


good. 
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FIGURE 18. Pressure — Temperature — Enthalpy Diagram for 42.3 per cent of Methane in Carbon Dioxide 
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FIGURE 19. Comparison of Excess Enthalpy for Methane — 
Carbon Dioxide Mixtures 
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FIGURE 20. Comparison of Excess Enthalpy for Methane 
— Carbon Dioxide Mixture 
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Experimental Results With Equimolal Methane - Carbon Dtoxide - 
Nitrogen Mixture 


The measurements of the isothermal throttling coefficient were 
made on an equimolal mixture of methane, carbon dioxide and nitrogen. The 
pressure range of the measurements was from 300 to 2000 psia and the three 
temperatures were 10, -10 and -30°C. The composition of the mixture was 
determined by the gas chromatograph and is given in Table IX. 

The smoothed isothermal data is tabulated at Table X and is 


Shown in Figure 21. Values of ¢ at zero pressure were calculated from 


the virial coefficients given by Dymond and mich an and me The 
calculation of $¢° is shown in Appendix C. 
TABLE IX 
Composition of Equimolal Methane - 
Carbon Dioxide - Nitrogen Mixture 
Composition Mole Fraction 
Me thane J2333 
Nitrogen Os353 
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Compartson_of Experimental Results 

Isothermal enthalpy departures for the mixture were calculated 
by integration of isothermal throttling coefficients with respect to 
pressure. The results were used to study the possibility of predicting 
ternary or multicomponent mixtures by using the binary experimental data. 
The experimental binary excess enthalpy data are from Lee and fee Oe 
for carbon dioxide-methane system and nitrogen-carbon dioxide system. The 
experimental excess enthalpy data of nitrogen-methane system are from Klein 
eepaied Pure component enthalpy data of nitrogen and methane are from 
Dice and carbon dioxide is from Vukalovich and paren 

Due to the limitation of the experimental excess enthalpy data 
available, the comparison is made only on isotherm of 10°C. The enthalpy 
departure of each equimolal binary system at 10°C and several different 


pressures was calculated from data of experimental excess enthalpy of the 


binary system and enthalpy of the pure component by the following relation, 


E 1 Uae ° 
= —_ -_ — 2 
5 =e aE 7 (H ae del) 5 (HE “teeiew) (28) 


where i is molal enthalpy departure of equimolal mixture of i and j 
he is molal excess enthalpy of equimolal mixture of i and j 
isle are molal enthalpy of i and j respectively 

and Hy > He are molal enthalpy at zero pressure of i and j respectively. 

The molal enthalpy departure of the equimolal ternary mixture can be 


calculated as follows: 


AH = =(AH,, + AH (29) 


HT,P 7 36 eHi2 + AHbs + His) pp 


The molal enthalpy departure of the equimolal ternary mixture of methane, 


carbon dioxide and nitrogen was also calculated by using BWR equation of 
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state with Bishnoi and Robinson mixing rules and parenerers 

The comparison between experimental results of this work and 
values calculated from binary data is shown in Figure 22. The prediction 
by BWR equation of state is seen to be in excellent agreement with these 
experimental results. The lack of agreement between this experimental 
data and the values calculated from binary data may be due to not taking 


account of ternary interactions. 


kapertmental Results With Equimolal Methane - Carbon Dioxide - 


Ethane Mixture 

Data of isothermal throttling coefficients were obtained at six 
isotherms for the equimolal mixture. The composition of the mixture 
which was determined by gas chromatograph is given in Table XI. The 
pressure range was from 250 to 2000 psia and the six isotherms are -10, 


0, 20°, 40, 60 and 90°C. 


TABLE XI 


Composition of Equimolal Methane - 


Carbon Dioxide - Ethane Mixture 


Composition Mole Fraction 
Methane 0.340 
Ethane ORB25 
Carbon dioxide O65 SSS 


The raw experimental data was treated by computer reduction 


method and the smooth isothermal throttling coefficients are tabulated in 
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Table XII. The smooth isothermal throttling coefficient at constant 
temperature is plotted as function of pressure and it is shown in Figure 
23. The zero pressure isothermal throttling coefficients were determined 


20,245 277%, 31,63,5102 


from second virial coefficient by using Equation (14). 


The calculation of $° is shown in Appendix C. 


Compartson of Results 


Peterson and Wilson reported isothermal enthalpy departure from 
20 to 2000 psia for this mixture at temperature from -50°F to 300°F. 
Enthalpy departures were calculated from the experimental isothermal 
throttling data which are listed in Table XII by integration with respect 
to pressure. The values obtained were plotted as function of temperature 
at constant pressure. Figure 24 shows the comparison between experimental 
data of Peterson and eis” and this work. The agreement is seen to be 


very good. 


Experimental Results wtth_50 Percent Methane in Carbon Dioxide 


and Ethane Mixture 
Isothermal throttling coefficient determinations of this mixture 
were obtained in the gaseous region at temperatures from -10 to 90°C, at 
pressure from 250 to 2000 psia. The composition of this mixture was 
determined by gas chromatograph and the compositions are given in Table 
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TABLE XII 
ISOTHERMAL THROTTLING COEFFICIENTS FOR THE 


EQUIMOLAL CO., = CH, CoHe MIXTURE 


4 
-¢ (J os fare) 


Temperature, °C 


PRESSURE =10.0 0 20 40 60 90 
(bar) 
Ox 41.80 40.20 34.60 30.40 26.80 22.60 
10 46.83 42.87 36.39 31639 IB) a5) 220 O2 
20 eos 47.40 38.99 32.96 28.46 23.41 
30 O32 54.16 42.18 34.74 MASS 23598 
38 17.29 - - - = 2 
40 - 64.55 46.22 36 .66 30.67 24.042 
50 = 80.96 1.50 39.01 31562 24.86 
60 - 114" 5 58.09 41.64 62.97 2D eae 
70 = 1422 65.04 44.16 34.06 ZN) 
80 = OOROL 70.61 46.13 34.91 25.83 
90 = 48.30 ORD D 47.17 35.234 25.90 
100 - 29.57 60.75 46.76 85.19 25h 
110 = Eas) 46.50 44.34 34.46 25.43 
120 - 14.30 34.92 3930 33.29 24.89 
130 = 10694 26925 34.36 Slee hey 2brelS 
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FIGURE 23. Isothermal Throttling Coefficient for the Equimolal Mixture 
of Methane, Carbon Dioxide and Ethane 
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FIGURE 24. Enthalpy Departures for the Equimolal Mixture of Methane, 
Carbon Dioxide and Ethane 
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TABLE XIII 
Composition of 50 Percent Methane in Carbon 


Dioxide and Ethane Mixture 


Composition Mole Fraction 
Methane 0.514 
Ethane 0) 6 Ae? 
Carbon Dioxide On244 


Isothermal throttling coefficients derived from the experimental 
data by computer reduction method are tabulated in fable XIV. 2Asplot of 
the isothermal throttling coefficient at constant pressure as a function 


of temperature is shown in Figure 25. 


Experimental Results wtth Equimoltal Methane - Carbon Dtoxtde - 
Nitrogen - Ethane Mixture 


Isothermal throttling coefficient determinations of this mixture 
were obtained in the gaseous region at temperatures from -30 to 60°C, at 
pressure from 250 to 2000 psia. The composition of this mixture was 
determined by the gas chromatograph and the compositions are given in 


Table XV. 
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TABLE XIV 
ISOTHERMAL THROTTLING COEFFICIENTS FOR THE 


50 PERCENT CH, IN CO, - C,H, MIXTURE 


4 1s 240 


= sell bene 


Temperature, °C 


-10 0) 20 40 
-20 32.80 28.80 25.20 
62 Sig 2 6s) 30.25 ZOe3 
83 Si cise) 31.67 Pier. 
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= O7, 66.39 42.60 Sve As 
~/3 10253 44.28 327.60 
2 64.61 44.89 S09 
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FIGURE 25. Isothermal Throttling Coefficient for the 50 Percent Mixture 


of Methane in Carbon Dioxide and Ethane 
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TABLE XV 
Composition of Equimolal Methane - Carbon 


Dioxide - Nitrogen - Ethane Mixture 


Composition Mole Fraction 
Methane 022535 
Carbon Dioxide Os255>D 
Ethane 0.2365 
Nitrogen 0.2545 


Smooth isothermal throttling coefficient of this mixture by 
this experiment are tabulated at Table XVI. A plot of the isothermal 
throttling coefficients at constant temperature as function of pressure 
is shown in Figure 26. Zero pressure isothermal throttling coefficients 


2092427, 31564570,101 


were calculated from second virial coefficients by 


Equation (14) and the calculation of $° is shown in Appendix C. 
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FIGURE 26. Isothermal Throttling Coefficient for the Equimolal Mixture 
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SECTION IV - EVALUATION AND EXTENSION OF METHODS OF PREDICTION 

This section presents a review of the methods of prediction of 
thermodynamic properties and particularly enthalpy of mixtures. Two new 
methods are proposed for the enthalpy prediction. 

To calculate the thermodynamic PEOperE1es Of aetiuld er fiuid 
mixture from a knowledge of molecular parameters by directly solving the 
equations of statistical mechanics, requires very detailed mathematical 
calculations as well as knowledge of the interactions between many 
particles. For this reason, no satisfactory solution has been achieved 
for all regions of enthalpy prediction. Statistical mechanics has been 
applied with some success in this endeavor especially with respect to the 
behavior of gases containing relatively simple molecules. 

Although advances have been made in the area of predicting 
macroscopic properties, it appears that for some time to come less 
sophisticated methods of prediction must play an important role. To obtain 
quantitative representation of thermodynamic properties, assumptions are 
required which make use of experimental data and certain parameters or 
functional forms having no direct fundamental significance. These 
approaches and other methods developed from correlating thermodynamic 


data will be discussed. 


Method of Prediction Based on Thermodynamte Data 


A recent review of the available methods of prediction of 
enthalpies of fluid mixture at elevated pressure was made by facnan: 
Mather?’ and verepES |p The main groups into which the methods can be 
divided are: 


1. Estimation of partial enthalpies 
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2. Equivalent pure component method 


3. Application of PVT data 
4. Application of equations of state 


5. Generalized corresponding states correlations 


L. Estimation of Parttal Enthalpies 

The enthalpy per mole of a mixture, H? can be determined 
exactly from a knowledge of the partial molal enthalpies of the individual 
components H, 5 by application of the expression 

H 


= z x Hy (30) 
il 


The partial molal enthalpy of a component is generally a function 
of composition. For gases at zero pressure, the assumption that the 
enthalpy of a component in a mixture is the same as the enthalpy of the 
pure component is an accurate one. Therefore Equation (30) can be applied 
rigorously to establish the enthalpy of gaseous mixtures at zero pressure. 
In general, however, partial molal data are not abundant, and therefore, 
several methods of estimation have been suggested. 

Very often the assumption is made that since hydrocarbons in a 
homologous series, are chemically similar, their solutions are ideal. 
Hence to compute a mixture properties, the sum of the property of each 
of the pure components (at the same temperature and total pressure) is 
taken in proportion to the quantity of the components present. This 
assumption may give good results at moderate pressures and relatively high 
temperatures. However, extreme care must be taken when applying this 
procedure since as Neenere has shown extreme values of the heat of 


mixing do exist, especially in critical regions. Charts for the pure 
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component enthalpies are given by Meenceubl pos Bapeeaae and Scheibel and 
mae A difficulty arises in using these charts because it is not 
always possible to obtain the enthalpy data for the pure component for the 
conditions of temperature and pressure under which the mixture exists. 
Peters has derived auxiliary curves for estimation of partial enthalpies 
under those conditions and Maxwell used an extension of the vapor pressure 
curve for the partial enthalpy of a low-boiling component in a liquid 
mixture. In general, the plots presented by Peters were limited to a 
temperature range between -260 to +420°F at pressure up to 600 psia. 
Similarily, Maxwell's plots extend between -200 and +200°F at pressure 
below 150 atm. 


A major drawback of this method is the fact that an enthalpy 


diagram must be available for every component present in a mixture. 


2. Equivalent Pure Component Method 


Several methods have been developed in which a new parameter is 
introduced to allow for the effect of composition. The data are usually 
presented as a series of enthalpy-temperature plots, there being a separ- 
ate one for each pressure, with lines of this constant parameter. The 
forms that this parameter takes are most commonly the molal average boiling 
point, the average molecular weight or the mixture specific gravity. 

Scheibel and enya presented monographs based on the average 
molecular weight of a hydrocarbon mixture. 

Papadopoulos etaalae showed that the molal average boiling point 
of a mixture defined as 


MeA.B.P. = 2 x, (B.P.) 5 (31) 
al 
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where the (B.P.), are the pure component boiling points, for lighter 
hydrocarbons served to correlate values of a partial molal enthalpies 
calculated from an equation of state. At the same time Canjar and 
Edmister’> employed the same correlation using the fugacity coefficients, 
employing a numerical technique for differentiation, and derived partial 
molal enthalpies for lighter hydrocarbons. Canjar and Papakea” prepared 
plots of the isothermal enthalpy departure as a function of temperature 
at pressure for mixtures with different molal average boiling points. 

The plots of Canjar and Peterka are restricted to the temperature range 
-200 to 500°F at pressures below 1500 psia for mixtures with molal 


average boiling points from -270 to 190°F. 


5. Application of P-Y-T Data 


Accurate P-Y-T data can be used to calculate the effect of 
pressure on enthalpy of fluids and fluid mixtures. If the volumetric 
data are available, the enthalpy can be calculated by use of the relation, 


OV 


P 
H(T,P) = H(T,0) + [. [V = T=), Jap, (32) 


where H(T,0O), the enthalpy at zero pressure, can be determined for most 
simple fluids. This is evaluated from ideal gas enthalpies of pure fluids 
either from statistical mechanics or from measured data. The ideal 
mixture enthalpy is determined by applying Equation (30). 

The term in brackets of Equation (32) under the integral sign 
involves the difference between two terms, one of which includes a 
derivative. As a result, extremely accurate volumetric data are required 
to yield reliable estimates of the effect of pressure on enthalpy. A 


reduction of accuracy of one order of magnitude is to be expected. 
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Volumetric data for mixtures of the required accuracy are available but 


very rare. 


4. Application of Equation of State 


A large number of equations of state have been proposed but only 
a few have been used for extensive calculations of enthalpies. In general, 
most equations of state serve to relate pressure P as the dependent variable 
to temperature, T, and specific volume V as independent variables. As a 
result, it is convenient to transform Equation (32) to, 


: P 


H(T,«) ~ H(T,V) = RT - PY + is (Pp - TC) Jay (33) 


Srv 20 


jee) 


The virial equation of state has a sound theoretical foundation and is 
free of arbitrary assumptions. The relations between virial coefficients 
and intermolecular potential functions can be derived from statistical 
mechanics and the composition dependence of the coefficients is known 
exactly. This equation is restricted to vapors since the power series 
relation tends to diverge at densities approaching those of liquids. 
; . , : ; V3 : 
Relatively simple equations like the Redlich-Kwong equation of 
State have been used for enthalpy predictions. Several modifications 
oar : 91 
had been attempted to improve the prediction, such as those by Wilson , 
76 F sy) i 
Soave and Chaudron, Asselineau and Renon  . More complex equations 
F Aa) : 56 
such as the Benedict-Webb-Rubin and Martin-Hou have also been used to 
predict thermodynamic properties, including enthalpy. New mixing rules 
? ale 
which were proposed by Bishnoi and Robinson for the BWR parameters are 
claimed to improve the accuracy of such predictions. The modification 


Dogo sel ah 


of BWR equation by Starling is claimed to improye the prediction 


especially at lower temperature regions. Other modifications of BWR type 


« ge 
ra a Nag ree, . ee 
vos 2 sits 0 ph aa Tia a 


laxdhes ot (aubytetind a0 baie . up gone: - 


aideivny FrahasGeb os aie ts eee 
sit A ee af 


s th _. detdbasee wn er : 
ov (SE) wnrieug? wana o3- > a a 


ae ma 7 ye 
as ait Er 9 ate bs tha = ss 


3} hit polishayot fastrotosata bevos # oof erage so suanips Intl 


_ 
_ 
a 


_ 


he  § aeer 


syretrdtieds Feisty nsdielsd eautielsx ad? “senoksquiees erat ieee 
caliebinge aT satin tk aay ping. Aaa 
imond at efustst3ia0n Spt Ye’ ayastimageb acksteagwen sd ee 7) 
eeittee xewog et sunita aungav 02 beaakizeet 21 mvizdeps eld yee 
Lisdup ti to sneid giAiasOaggh sotatndeb Ue amreetd oa saoliaas 
te «otieupe ©" gaowi-fob lind SAd svt nvolanips slants ylovigelet oe 
enstsesttibom)taravat .ehotiakseag yqindsas 402 banu amas und SIR 
soar td seqtia ex 4ibua jnoksokbeeq als Svorqat 07 bstqmsoa asad bad 

aootsiups-atqncs, tom. iwnal beh cedatioeas’ , ocviubity Bao” meee 

pa haa ated oath ovat “varanasi Bae “abtint-ddweaslienee a9 ee Aue 
esfut gniate wat Jghadghe, guipulanl ,absrsacns sinien(bemtedd 22tbe7¢.. 

sxe, Hesonnaug ideas 403 apanildae bie. tondethl vd islingasg wade a5 
fotsext than of ieheticthing dove Ye yasrose adv aevaged Os Benita 
nakvabiionar ots ovovgar ox baat at 8" oc ragaet 9B weteyps MME Tos 
sige? mult tormaei sant? Uhl reds) aridity ~ancitags « Silliae 


81 
Of equations of “state f"such as@that pede, and Mecareyae have also been 
proposed. Recently parameters for BWRURIe ieee mi ak have been 
determined by multiproperties regressing on volumetric data, vapor 
pressure data, enthalpy and k-value data. This ensures thermodynamic 
consistency of all type of data. 
Unfortunately neither accurate thermodynamic data nor 


parameters for equation of state such as BWR which need large amounts of 


experimental data are available for many components. 


0. Generaltzed Corresponding States Correlations 

Recently reviews of the application of the corresponding states 
principle have been presented by setae and Leland and Ghtone1eane 
These methods are based on the theory of corresponding states first put 
forward by van der eee in which he stated that all pure gases would 
have the same compressibility factor when measured at the same (reduced 
condition of) temperature and pressure. The law of corresponding states 
can be stated mathematically 

De ie gs een Rama (34) 

J i J A 
where the symbols Yi» and vs represent corresponding dimensionless group 
Went Wan 


properties of two fluids "i'' and "j" given as function of their respective 


reduced temperature aR 3 ies and reduced pressure ae = ace 


This statement of principle has become known as Two-Parameter 
Law of corresponding state because one must have knowledge of the two 


parameters (T,, and Po) in order to apply it and must assume that there are 


C 


only two adjustable parameters in the intermolecular potential in order 


to derive it from statistical mechanics. 
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Lydersen, Greenkorn and Hougen (LcH)?4 correlated the 
thermodynamic properties of pure compounds based on a modified correspond- 
ing states principle. This principle states that pure compound with the 
same critical compressibility Zo values, the dimensionless groups of 
properties of two pure compounds have the same numerical values at 
corresponding reduced condition. The correlation of LGH has the form 


H(T,0) - H(T,P) Ho - 4H 


= [———] + (Z., - 0.27) [D] (35) 
i as C 


where the bracketed terms are presented as generalized functions in 
tabular form by the authors. The condition covered include Pa =< 305for 

95 ; 96° 
0.5 <T, < 15. Yen and Alexander ~ and Yen and Garcia-Rangel” ~ improved 
the correlation of the isothermal effects on enthalpies for pure compounds 
with the aid of reliable literature data. 
PreZen die ae employed the acentric factor w, which is related 


to the shape of the reduced vapor pressure curve as the third parameter. 


1 : 
The correlation of Curl and Pitzer ; is given in the forn, 


0 Cy) 


H(T,0) ~ H(T,P) fie ate (eed 4 
(36) 
RT, RT RT, 


where the bracketed terms are somewhat different generalized function 
presented in tabular form by the authors. These tables cover the range 
of pressure for PR Se nape W's) <T, < 4. Revision of the original 
correlations which incorporated enthalpy data at elevated pressure in 


94 
addition to PVT data have been presented recently (Yarborough ). Chao 


and eyseplysea = recently extended the Curl and Pitzer correlation tempera- 
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ture range down to TR = 0.35. The extended enthalpy correlation was 
obtained upon differentiation of the generalized fugacity function 
previously developed by Chao et al?® and upon combining with calorimetric 
data in the literature. 

Leach and Wina e extended the two parameter law of 
corresponding dates to include fluids which have structural dissimilarity. 
Molecular shape factors were introduced which modified the critical 
temperature and critical volume of fluids whose molecular structures were 


dissimilar to those of a fluid for which thermodynamic data were 


available. 


Molecular shape factors (85, and Soe) Of a fluid i relative to 


fluid "j" are defined mathematically by the equations 


Thy, Ihe (37) 
~ i 
AH AH 
a) ee ec (38) 
RT, a RT, 5 


at corresponding modified reduced conditions, 
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where Z is the compressibility factor of the two fluids. 
In general, the molecular shape factors of a fluid are functions 
of its temperature, density and intermolecular parameters. Leach 


correlated the shape factor for isomeric and paraffin hydrocarbon using 
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methane as a reference. Watson and Row brisent: and Gunning and 

; 29 F 
Rowlinson used the same correlation for non-hydrocarbon systems. The 
molecular shape factor can be extended to predict thermodynamic properties 


of the mixtures. The mixing rules have been suggested by tena”. 


Appltecatton of Generalized Corresponding States Correlattons to Mixtures 


In applying the above correlation to mixtures, it is assumed 
that the reduced functions for a mixture behave in the same manner as 
those for a pure component. It is necessary to establish values of the 


three parameters for mixture. The concept of pseudocritical properties 


was first suggested by ea The suggested linear form are: 
cone: dat cs ce 
my i i 
ey eae (42) 
m sl af 
and 
Oe ——uax Oy, (43) 
m , ii 


There are several other mixing rules which have been proposed in the 
literature. The procedure and the original background of the derivation 


of the mixing rules is presented. 


69 
4) Method of Pitzer and Hultgren 


Pitzer and Hultgren had extended the three parameter acentric 
factor corresponding state theory to binary gas mixtures, using the. 


pseudocritical concept. By studying the experimental volumetric data, 
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it was found that simple quadratic formulas could be used to represent 
the composition dependence of the three pseudoproperties of each mixture. 


Applied to the multicomponent mixture, the equations of the mixing rule 


are: 
Le ee ee oe) 
m ie} ale 
P cligited Xj ique (49) 
m cat ley 
Shen xe x0 (46) 
m lop) mreigie Pasi | 
1 J 
The constants To 5 Po and Oe (i#}), must be evaluated from data on 
14 ey 


binary mixtures. These mixing rules reduces to Kay's rules when the 

constants To : Po and W,, are taken as arithmetric average of the pure 
ij ij J 

constants. 


it) Method of Joffe-Stewart, Burkhardt_and (OO ee 


The equations of the mixing rules are given as, 
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in van der Waals equation where, 


an i" " Yi Oe) 
1 

b = 

m F ; ¥4¥ Pay (52) 


Since the van der Waals equation is a two parameter equation, the strict 
development is limited to mutually conformal substances. However, the 
rules can be applied to other mixtures by assuming a linear variation in 


the third parameter. 


tit) Method of Leland_and Meter 


The problem of defining pseudocritical values for a multicompo- 
nent phase is analogous to the problem of writing the virial coefficients 
for the pure components. The second virial coefficient for an m-component 


mixture is related to the pure component virials by: 


Bees ey iyes (53) 


Ae eee ij 

ij j ij 
Leland and Mueller obtained the mixing rules by equating terms of a second 
virial coefficient expression for mixtures with those of a pure component 
and making appropriate simplifying assumptions. The pseudocritical 


constants are given as 
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(55) 


and is tabulated by the authors. The pseudocritical constant obtained by 
this method is not only function of components of mixture but also 
function of temperature and pressure. Again the third, parameter is 


assumed to be linear variation for the mixture. 


tv) Method of Prausnitz and Pot ee 


The mixing rules of Prausnitz and Gunn are also suggested from 
relations for the second virial coefficient of conformal substances. 


These rules are given as: 
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(59) 


(60) 


(61) 


(62) 


(63) 


(64) 


where the A terms are small correction terms depending on the nature of 


the binary pairs. 


The quantities r and s are functions of TR 


tabulated by the authors. 
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Prausnitz and Gunn also recommended the following simplified 
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RT 
Cn ay a ce (65) 
m f, i C; i all 


where again Za is given by Equation (60). 
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v) Method of Barner and Cuan 

Barner and Quinlan suggested the mixing rules based on a modi- 
fied form of Redlich-Kwong equation of state (Chueh and Prsdenite or The 
mixing rules for the two constants of the equation of state is given by 


Joffe and Badkewittchine as, 
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The suggested pseudocritical constant of the mixture is given by 
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The interaction constants sig are tabulated by the authors. coe have been 
generalized as a function of Vo [Vo and are also given by the authors. 


Le 


vt) Modtfted Method of Gunn 


enties developed a corresponding states theory for fluid mixture. 
The pseudocritical temperature is obtained by using generalized equation 
for second virial coefficient of Curl and Sines © and the mixing rule for 


second virial coefficient. 
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The rigorous expression for the scaling temperature is obtained by equating 


Pauation. (53) and (71). 
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A new correlation of second virial coefficient of both polar 


and non-polar system presented by aaeneocr ioe te replaced that of Curl and 


: 70 ; ; 4 
Pitzer . The new correlation for non-polar system is given by, 
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The interaction constant Si is given by Prausnitz or by fitting the 


C 


experimental Boa with Equation (74). The pseudocritical constant T, and 
m 


Po is obtained by solving Equation (77) and (68) simultaneously. 
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vit) Proposed Method (Soave-Redltch-kwong Equation of State) 


A modified Redlich-Kwong equation of state has been proposed by 
see. The parameter a in the original equation is suggested to be 


temperature dependent. The modified equation is, 
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Applying the equation of state to the mixture, the mixing rules for the 


two parameters are, 
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where ie is an empirical correlation factor to be determined from the 

experimental vapor-liquid equilibrium data, for each binary component. 

present in the mixture. Generalization of Se has been arrempredon 
Starting by substituting Equation (82), (84) and (88) into 


Equation (66) and (87) respectively, 
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To and Po are obtained by solving Equation (91) and (92) simultaneously. 
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Here the linear variation for wo is assumed. 
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Evaluation of Nine Enthalpy Correlations 

The nine correlations evaluated and compared in this study are 
shown in Table XVII. Six of the correlations used corresponding states 
method of Curl and Bieeer with six different mixing rules and three of 
the correlations used equations of state. The following six mixing rules 
for mixtures which have been discussed above to estimate the pseudo- 
critical parameters for the Curl and Pitzer generalized correlation are: 
1) method of Barner and Quinlan, 2) proposed method (Soave-Redlich- 
Kwong equation of state), 3) modified method of Gunn, 4) method of 
Prausnitz and Gunn, 5) method of Joffe-Stewart, Burkhardt and Voo, and 
6) method of Leland and Mueller. The equation of state correlations are 


jaiad . SoavesRed Prene kwon. - equation and Mark wee correlation. 


Starling-BWR 
; : 1, ; 

In this study, the Curl and Pitzer generalized correlation 

tables for enthalpy departure and the extended correlation temperature 
range of Chao and Creenicarnas were used. The linear interpolation, 


performed by computer was utilized in the calculation of the enthalpy 


departure. 


Presented in Table XVIII are the data sources, number of data 
point, and temperature and pressure regions of the data for the sixteen 
fluid compositions studied in this evaluation. 

Presented in Table XIX for each correlation are the average 
absolute deviations between predicted enthalpy departures and experimental 
enthalpy departures for each of the sixteen fluid compositions used in 
this study. The overall average absolute deviation for the total data 
points is given for each correlation in Table XIX. It is to be noted 


that the BWRS yields the lowest average absolute deviations. Results 
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TABLE XVII 


CORRELATIONS EVALUATED IN THIS STUDY 


Correlation 


Starling-BWR 
Soave-Redlich-Kwong 
Mark V 

Corresponding State + 


Method of Barner & 


Quinlan 
Proposed Method 


Method of Prausnitz & 


Gunn 


Method of Joff-Stewart, 
Burkhardt & Voo 


Method of Leland & 
Mueller 


Modified Method of 


Gunn 


Abbreviation 


CPBQ 


CPNM 


CPPG 


CPSV 


CPLM 


CPMG 


Developers 


Starling 
Soave 
Wilson 

Curl, Pitzer 


Barner, Quinlan 


This work 


Prausnitz, Gunn 
Joffe,Stewart, 
Burkhardt, Voo 


Leland, Mueller 


Gunn 
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obtained using method two listed above, gives the lowest average absolute 


deviations among the six methods. 


The pseudocritical parameters (temperature and pressure) which 
were obtained by the five different methods for the mixtures of 49.4 mole 
percent methane 7 propane plot as function of temperature are shown in 
Figure 27 and 28. Pseudocritical temperature and pressure obtained by 
modified method of Gunn shows very large deviation from others. Pseudo- 
critical temperature and pressure obtained by method of Leland and 
Mueller are the function of temperature and pressure. They can not be 


shown in Figure 27 and 28. 


Conelustons 

The conclusion of this study is the finding that the equations 
of state BWRS, Soave-RK, Mark V and the corresponding states Wieh tire 
proposed mixing rule which is derived from Soave-RK, are the most 
accurate of the nine enthalpy correlations tested. Conclusions relating 


to each correlation individually are given below. 


Starling-BWR 

The average absolute deviation for this correlation is 0.9 Btu/ 
lb. This method is not a generalized correlation and therefore can be 
applied only when the parameters of this equation of state have been 
determined for each component in a mixture. Also binary data are needed 


to determine interaction constants. 


Soave-Redlt ch- kwong Equation of State 


The average absolute deviation for this correlation is ean 
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Btu/lb. This method can be considered as a generalized correlation because 
only critical properties and acentric factor of each component in mixture 
are required for the prediction. The only additional information is the 


interaction constant k, which has to be determined from the binary 


j 


vapor-liquid equilibrium data. Generalization of bie for some system 


6 : : 
has been attempted o This method requires very little computation time. 


Mark V Computer Program 
The average absolute deviation for this correlation is 1.60 
Btu/lb. This correlation always predicts too high at low pressure and 


predicts too low at high pressure. 


Corresponding State of Curl and Pitzer With Different Mixing Rules 


Proposed Method 


The average absolute deviation for this correlation is 1.70 
Btu/lb. The interaction constant 5 for the Soave-Redlich-Kwong equation 
of state has to be determined from the binary vapor-liquid equilibrium 


data. The generalization of oe, has been attempted. 


Method of Barner and Quinlan 

The average absolute deviation for this correlation is 2.0 Btu/ 
lb. The mixing rules is derived from modified Redlich-Kwong equation of 
state (Chueh and Peausnitee Barner and Quinlan had generalized the 
binary interaction constant Ra for most of the system. The pseudocritical 


temperature and pressure of the mixture are independent of the temperature 


of the system. 
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Method of Joffe, Stewart, Burkhardt and Voo 

This correlation gives average absolute deviation of 2.32 Btu/ 
lb. The mixing rules are derived based on van der Waals equation of 
state. No binary interaction constant is needed. It is not surprising 


that the prediction is less accurate than other methods which have one 


more adjustable parameter. 


Method of Leland and Mueller 

The average absolute deviation for this correlation is 2.6 Btu/ 
lb.. The pseudocritical temperature and pressure obtained by this method 
are both functions of the temperature and pressure of the system. This 
correlation gives very good prediction for some systems but gives very 


poor prediction for some other systems. 


Modified Method of Gunn 

This correlation gives very poor prediction in liquid region, 
but gives good prediction in gas region. The experimental data of this 
work are in gas region. It is shown in Table XIX that this correlation 
gives the lowest average absolute deviation among the corresponding states 


methods. 


Method of Prausnitz and Gunn 


The average absolute deviation for this correlation is 1.91 


Btu/1lb. Equation of the form Bo = 22 eye forms the theoretical 
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basis on which the prediction of pseudocritical constants is based. The 


two empirical constants ATo and AV, for each binary had been 
ij 1j 
generalized by the authors. Two more empirical constants s andr for 


each mixture are needed to obtain the pseudocritical constants for the 


mixture. 
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Summary and Conelustons 


An isothermal throttling calorimeter for direct determination of 
the effect of pressure on enthalpy was designed, fabricated, and 
tested. The recycle flow facility, as described by Bishnoi, was used 
to recycle the material for this study. 

The effect of pressure on enthalpy was determined experimentally 
for nitrogen at temperatures from 0 to 100°C and pressures from 300 
psia to 2000 psia. The data are in good agreement (within + 1%) with 
other direct determinations in recent literature. 

Measurement of the effect of pressure on enthalpy were made on 
14.5 mole percent and 42.3 mole percent of methane in carbon dioxide 
at gas region at pressure from 300 to 2000 psia. 

Skeleton enthalpy tables and enthalpy-pressure-temperature 
diagrams for those two mixtures were prepared, using isothermal 
enthalpy data of this work and isobaric heat capacity data of Bishnoi. 

Measurements of the effect of pressure on enthalpy were also made 
on equimolal mixture of methane, carbon dioxide and nitrogen, 
equimolal mixture of methane, carbon dioxide and ethane, methane rich 
mixture of methane, carbon dioxide and ethane, and equimolal mixture 
of methane, carbon dioxide, nitrogen and ethane. The measurements are 
in gas region and pressure up to 2000 psia. 

The data obtained in the course of this investigation plus data in 
the literature were used to compare several of the available methods 
of prediction. This comparative study indicates that corresponding 
states principle would be a most fruitful approach for extending 


methods of prediction. 


A new recipe for obtaining pseudocritical parameters for the 
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mixture of the three parameter corresponding states correlation of 
Curl and Pitzer was developed to predict the isothermal enthalpy 


departure of mixture. In comparison with sixteen fluid compositions 


of experimental data in the literature, this new method gives lowest 


average absolute deviation among other methods of obtaining pseudocri- 


tical parameters. 
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Recommendations for Future Work 


The pressure transducer for measuring the pressure drop across the 
calorimeter should be replaced by a more accurate and wider range 
pressure transducer. 

The new device of flow meter should be developed in order to 
obtain more accurate experimental data. 

The calorimeter should be modified in such a way that it would 
be easier to interchange the throttling coil. 

A higher pressure drop throttling coil should be constructed in 
order to study the effect of pressure on enthalpy in liquid region 
and region across the two phases. 

An oil and water removal unit should be installed at the outlet 
of the compressor in order to prevent the solid formation of o1lF or 
water which will plug the throttling unit of the calorimeter when 


operated at low temperature. 
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NOMENCLATURE 


Constant in Equation (21) 

Constant in BWR or BWRS equation of state 
Constant in RK or SRK equation of state 
Constant in BWR or BWRS equation of state 
Calibration constants (flow meter) 

Second virial coefficient in volume expansion 
Second virial coefficient in pressure expansion 
Reduced second virial coefficient 

Normal boiling point 

Third virial coefficient in volume expansion 
Third virial coefficient in pressure expansion 
Isobaric heat capacity 


Correlation term of Lydersen, Greenkorn, Hougen 


correlation 
Base of natural logarithm 
Mass flow rate 
Specific enthalpy 
Partial molal enthalpy 


Expression in mixing rule of Joffe-Stewart, Burkhardt, 


Voo 


Interaction constant for Barner and Quinlan mixing 


rule 
Interaction constant 
Molecular weight 
Slope of so against Tee Equation (85) 
Pressure 


Critical pressure 
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Reduced pressure P/P. 


Pressure drop 

Rate of heat transfer 

Work per unit mass W/F 

Gas constant 

Correction factor in mixing rule of Prausnitz and Gunn 
Temperature 

Critical temperature 

Reduced temperature T/T, 

Specific volume 

Rate of transfer of work 

Mole fraction 

Mole fraction 

Compressibility factor 

Constant in mixing rule of Leland and Mueller 
Constant in BWR or BWRS equation of state 
Expression in mixing rule of Prausnitz and Gunn 
Difference 

Joule-Thomson coefficient 

Viscosity 

Density 

Dimensionless group of a thermodynamic property 
Molecular shape factors 

Reciprocal of temperature BA 

Acentric factor 

Parameter in RK or SRK equation of state 


Isothermal throttling coefficient 
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Subsertpts 

ad Adiabatic condition 

C Critical point property 

HAG Component in a mixture 

m Mean value 

m,mix Mixture property 

fe) Zero pressure value 

0 Reference substance property 
L Component 1 

ie Inlet condition 

uae Component 1 

ifs Interaction between component 1 and component 2 
B22 Component 2 

2 Outlet condition 


Superseripts 


Zero pressure value 


E Excess property 
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TABLE XX 


CALIBRATION DATA FOR 8163 MODEL PLATINUM RESISTANCE THERMOMETER 


This thermometer was calibrated for use in the range 
90.188 K < T < 773.15 K (-182.962°C < t < 500.60°C) on the International 
practical temperature scale of 1968. 


FOr gles .6 ly ke eee oak 
W = Wect + dw 


where W = R/R(0°C) and Wect is a reference function, independent of any 
particular thermometer. 


For. 90.186 Kk <9 1-2/5. 160k 
dw = A4 t + C4 (t — 100) eS 
For 0 Guc to 650.7200 
f= t edt 
where t' is determined by 
v 1x2 
w= 1+ At' + B(t') 
and 


i Lar Eee 
dt" = 0.45(t'/100) Gog - 1) Giovsa - 1) Gao77a - 


The coefficient in the above relations were found by fixed-point 
calibrations, i.e. triple point of water, the steam of tin point and the 


zinc point using continuous thermometer currents of 1 and 2 mA. 


yas to Insbasqsbnt om mn oe ro a> 7 


ici ee of 


. 


a a 
S°AU.0Ed > 2 > OO roe. 
Sa, 
ab+*2e 7 
> os 


vd Bankeyeseb ab "3 oxofiw 


S(aya + a+ bw 


A4 
C4 


R(0°C) 


i] 


TABLE XX (continued) 
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TABLE XXTI 


CALIBRATION OF HIGH PRESSURE GAUGE 
Manufacturer: Heise 


Serial No.: H39686 


Actual Pressure Heise Gauge 

(psi) oe down 

0 0 0 
200 200 200 
400 400 400 
600 596 600 
800 796 800 
1000 996 1000 
1200 1200 1200 
1400 1400 1400 
1600 1600 1600 
1800 1800 1806 
2000 2000 2005 
2200 2200 2205 
2400 2400 2405 
2600 2600 2605 
2800 2800 2804 
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TABLE XXII 


CALIBRATION OF 250 PSI DIFFERENTIAL PRESSURE TRANSDUCER 


Manufacturer: Validyne 


psi Voltage 
2uDe down 
0 0.001 0.004 
30 69 sopesne 
50 1.956 13976 
70 Za 2.760 
80 Jl 3D Slow 
100 Eye PAS) Bio: 
120 4.717 4.743 
140 5.514 5.541 
160 6.316 Ones) 
180 dilee 7.144 
200 15930 12952 
220 SaioZ 8.766 
240 Deoi a 902 


250 On 992 Teooe 
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CALIBRATION OF 75 PSI DIFFERENTIAL PRESSURE TRANSDUCER 


TABLE XXIII 


Manufacturer: 


Validyne 


Voltage (Volt 


0.001 
0.654 
Egea bil 
bese cae 
2.635 
3.303 
3916 
4.656 
4.999 
4.657 
3.980 
3.310 
2.642 
Lee) 
i316 
0.658 


0.001 
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CALIBRATION OF 60 INS, OF WATER DIFFERENTIAL 


TABLE XXIV 


PRESSURE TRANSDUCER 


Manufacturer: 


Inches of Water 


14 


Ly 


2293 
945 
»541 
“hae 
- 209 
O41 
-/0 
70 
cor 
~45 
SUD 
a2) 
. 60 
See: 
Beh) 
.043 
354 


aLOU 


310 


ole 


Foxboro 


Voltage (mV) 
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TABLE XXV 


FLOW METER CALIBRATION 


(Orifice) 
EQUATION: = Aces 
Leech 
where 128 gm-mol fa 
Agee br Calibration Constant 


Nee 
k CEG 
P 

y: = 227 2 

B: Ratio of orifice I.D. to tube I.D. 
NEYO Pressure drop across the orifice 
PB: Inlet-pressure of the orifice 

Ze Compressibility factor of the gas 


(gas mixture) 


P 
psia 


44.0600 
44.1100 
44.1100 
44.1100 
44.0100 
44.1100 
44.0100 
44.1400 
44.1400 
44.1900 
44.1600 
44.1600 
44.1400 
44.0 100 
43.8700 
43.8700 
43.8700 
43.8700 
43.8700 
43.8700 
43.8700 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 
43.9100 


AP 
psi 


1.8476 
1.5086 
Nee 35 
1.1880 
0.7536 
1.0314 
0.3266 
elitr 
Lute Movil lis) 
Peco 
1.0469 
Oem (sks: 
0.4266 
0.2342 
0.2860 
0.8162 
0.4917 
0.6184 
0.7413 
0.9240 
1.0988 
1.5861 
1.7790 
1.5407 
0.9288 
0.4239 
0.4426 
0.6237 
0.8364 
1.0149 
Lo? 1 
UI eheyexs) 
1.5428 
165634 
©. 10n2 
0.1343 
0.2459 
Oe 909 
O% 3389 
OMI 
0.2347 
0.2492 
Oem ets) 
0.4111 
262406 


TABLE XXY 


F 


4.8095 
4.3566 
4.0145 
3.8783 
3.0989 
366234 
2.0726 
Seusas 
424199 
4.0049 
326544 
21370 
2.3446 
e528 
202932 
226288 
Ze 9907 
3263540 
3.6774 
4.0825 
4.4605 
5.3041 
5.25864 
522467 
4.0697 
2.4950 
22-8506 
eRe cle ke 
Beet 
4.3044 
4.6666 
4.9209 
53 26 
Sip Shisie i) 
1.3637 
WA ree 
2 LGA 9 
1.8980 
Zeon.d2 
EwOONS 
22-0941 
lei229 
De Joo 
(A CAUSES) 
5.0162 
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= (ro oe 
-@.0130 
-0.0160 
-0.0167 
=O 30 128 
0.0107 
=) S0230 
-0.0242 
=) s00 7S 
-0 00104 
-0.0044 
-0.0124 
0.0053 
0.0046 
-0.0001 
0.0016 
0.0031 
0.0106 
-0.00 24 
0.0154 
-0.0068 
-0.0246 
0.0093 
O02 
0.0109 
0.0062 
0.0241 
0.0373 
0.0260 
0.0303 
-0.0034 
0.0308 
0.0468 
0.0000 
0.0141 
Os0 225 
0.0253 
0.0261 
0.0186 
0.0183 
0.0062 
0.0098 
0.0017 
=) ,0027 
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TABLE XXV (continued) 
P AP 
F ; YiPsAP 

psia psi 7M F Fal FOF al 
43.9700 22460 1.4490 5.0162 5.0247 -0.0085 
43.9700 220846 1.3991 4.8654 4.8517 0.0136 
43.9700 1.8618 1.3264 4.5785 4.5993 -0.0208 
43.9700 1.6395 1.2486 4.3140 4.3291 = .0b5 1 
43.9700 1.4188 be loSt 4.0276 4.0393 -0.0117 
43.9700 1.2644 Teo 3.8065 Soe wk -0 .0146 
43.9700 1.0811 1.0217 3.5360 325420 -0.0060 
43.9700 0.9208 0.9450 8.2599 S.2108 =(7..0)) 5 
43.9700 0.7589 0.8598 229649 2.9801 =O 152 
43.9700 0.6184 O.7776 22-6850 226948 —0,0098 
43.9700 0.5174 Oni 224418 224681 —0.0263 
43.9700 0.24453 0.6614 202184 262915 -0.0131 
43.9100 0.4490 0.6637 262055 222994 0.0060 
43.9100 0.6750 0.8113 22 1892 228116 -0.0224 
43.9100 0.8610 0.9139 3.1566 3.1677 -~0.0111 
43.9100 1.0768 1.0190 Benso0 3.532 0.0032 
43.9100 LAT MA) 1.1070 3.8289 3.8378 -0.0089 
43.9100 1.6475 1.2506 4.3140 4.3361 -0.0221 
43.9100 Ve9215 1.3474 4.6891 4.6721 0.0169 
43.9100 22-1274 1.4115 4.9146 4.8947 0.0198 
44.2400 0.3945 0. 1602 2¢ 6380 2 06345 0.0034 
43.7800 1.7865 1.5696 524560 5204433 0.0126 
43.7800 1.3953 1. 3950 4.8305 4.8373 -0.0068 
43.7800 0.08893 Lelie 3.8984 328889 0.0094 
43.7800 0.5265 0.8675 229951 3.0067 -0.0116 
43./800 0.3854 0.7437 225586 we ON Wd -0.0185 
43.7800 0.22433 0.5920 220415 2.0508 -0.0093 
43.7800 0.1663 0.4901 1.6862 1.6970 -0.0108 
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TOTAL SUM OF SQUARE 0.97146E 03 


SUM OF RESIDUAL = O.2823fE-05 


SUM OF RESIDUAL SQUARE = 9.20049E-01 


PERCENT STANDARD DEVIATION = 0.48473 


STANDARD DEVIATION = 9.16804E-01 


Y= A+ B¥X 
A= -0 .00 391 
B= 3470 33 
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Caltbration of Gas Chromatograph 


Ll) Methane - Carbon Dioxide Mixtures 

The binary mixtures of carbon dioxide and methane were analysed 
on a Burrell K-2 Kromo-tog Chromatograph using a three foot column of 
one-quarter inch diameter and packed with silica gel (40-100 mesh). A 
thermal conductivity detector was employed and helium was used as a 


carrier gas. Other variables were maintained at the following conditions: 


Carrier gas flow rate - 50 ml/min 
Column temperature - 80°F 
Detector Current - 150 mA 
Sample size = PL julit 


Excellent separations, with sharp well-defined peaks were obtained and 
the area under a peak was determined by a disc integrator. Standard 
sample mixtures of carbon dioxide and mixture were prepared by weighing 
small sample bottles and the amounts of methane and carbon dioxide 
charged. The weights of the empty bottle and the mixtures sample 
charged were approximately 155 g, and 1.75 g, respectively. The weighings 
could be reproduced to within + 0.1 mg. The error in a mixture composi- 
tion thus prepared should not exceed + 0.2 percent. The standards were 
then analysed on the chromatograph and the area fraction of methane (area 
under a methane peak divided by the sum of the area under the carbon 
dioxide and methane peaks) was correlated to the known mole fraction of 


methane in the standards by the following functional relationship: 


y = AK + BX 


where A, B are constants 
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tas 
Il 


area fraction of methane 


K 
i] 


mole fraction of methane 
The constants, A and B, in the above relationship were determined 


by a least square fit and were found to be as follows: 


> 
i} 


0.73451 


ee] 
I 


20 2720 


The percent standard deviation of the fit was found to be + G.47. 7A 


plot of Y versus X is shown in Figure 30. 


2. Methane, Carbon Dioxide and Ethane Mixture 

The mixtures were analysed on the same chromatograph using a 
one meter column of one-quarter inch diameter packed with porapak Type Q 
(80-100 mesh). The gas chromatograph was operated under the following 


conditions: 


Carrier Gas - Helium 
Carrier gas flowrate - 40 ml/min 
Column temperature - 80°F 
Detector current - 170 mA 
Sample size - 2ml 


Standard mixtures of methane — carbon dioxide - ethane were 
prepared in the same method as for binary mixtures. The standards were 
then analysed on the gas chromatograph and the area fraction of each 
component were plotted against its more fraction as shown in Figure 31. 
A straight line was drawn through each component. The response factors 


of each component were then determined from the calibration curve by 
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Mole Fraction of Methane 


Cas Chromatograph Calibration for Methane-Carbon Dioxide Mixture 
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TABLE XXVI 


RESPONSE FACTOR CALCULATED FROM CALIBRATION 


Methane - Carbon Dioxide - Ethane System: 


CH, : 36 
co, : 49 
CoH. : 53.5 


CH, : 36 
Ny 41.8 
co, : 49 


Methane - Carbon Dioxide - Nitrogen - Ethane System: 


CH, : 36 
co, 48.6 
CoH 54.0 


(Nexis) 
assigning the response factor of methane equal 36. The calibrated 


response factors were listed in Table XXVI. 


3. Methane - Carbon Dtoxtde - Nitrogen Mixture 

The same gas chromatograph with same conditions and same 
procedure of the previous ternary system were used to calibrate mixture 
of this system. Figure 32 shows the area fractions of each component 
plotted as function of its mole fraction, and the calibrated response 


factors are listed in Table XXVI. 


4. Methane - Carbon Dioxide - Nitrogen - Ethane Mixtures 

The same gas chromatograph with same operating conditions and 
same procedure of the previous ternary system were used to calibrate 
mixture of methane - carbon dioxide - nitrogen - ethane system. Figure 
33 shows the calibration curve and the response factors are listed in 


Table XXVI 
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Area Fraction 


Mole Fraction 


Figure 31: Gas Chromatograph Calibration for CH, -CO.-CoHe Mixture 
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Error Analysts 
The desired values of the isothermal throttling coefficient is 


determined from the results of measurement: 


OB, w sig oWe 
va = oles T°” AP F (Al) 


where W = power supply 


F 


mass flow rate 


AP pressure drop across the throttling unit 


The maximum relative error of Die in the experiment is 


Ag _ AW , A(AP) ag. 
o WenhPae eee ate PT x? m 


(A2) 


(oy 
Gon fy Co oi 


where AT error in the measurement of inlet temperature of the gas in 


the calorimeter 


AP = error in the measurement of the system pressure Pa 
Ax = error in the measurement of composition in the mixture 
Clearly, 
am DAW W 
AW H AW loss 
=e Sa — oe A 
ww, * CWrer 7H Ce) 


calculated error of measurement of the compensating power 


where AW, / Wy 
of the calorimeter heater; (AW/W) 37 = error in the value of the compen- 
sating power, due to the non-zero value of the temperature difference 
between the inlet and outlet; AW, |. ./W < error connected with spurious 


heat transfer between the calorimeter and the thermostat. 
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In this experiment, a Hewlett-Packard Model 3450 multifunction 
meter is used for all electrical measurements. For W = IV, the measure- 
ment of current I is by means of measuring the voltage drop across a 
standard resistor, the uncertainty of voltage is about + 5 microvolt. 

So the term (AW/W) 5 = +(0.005 - 0.002) %. The term (AW/W) 07 can be 
replaced by the ratio AT, / AT. 4> where AT, = the degree of non-isothermal 
state of the inlet and outlet of the gas when the calorimeter is in 
operation; AT od = degree of non-isothermal state when no heat is added. 
The deviation of null meter which is used for the measurement of the 
output of the five junction thermopile is about + 2 uV which is equivalent 


to O7017C) then 


ayer OS0l eee emery 

AT Wd AT Wd AT od 
For nitrogen: Out. < CAW/W) oy a ysis 
For CH, (43.3%)-CO, : O.03%> < CAW/W) a7 <OLe 


For Equimolal CO, CH, -CoH, O;03704 CAW/W) a7 Dally 
The term CAA is minimized by using the radiation shield surrounding 


the throttling unit and the calorimeter is evacuated to about 5 microns. 


This term can be neglected. The error in determining AP and F are: 
for NCAP) Se 40.5% 
AP ane 


and for (AF/F): 


(AF/F) < +0.82% for nitrogen 


< +1.02% for 42.32 CH, in co, mixture 


OM a7 for other mixtures 
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a2 
so, the maximum systematic error without the related error i.e. without 


allowing for the last term in Equation (A2) is 


ONO) 5S toy ta ualliee for nitrogen 
lO Fees ae ue for 42.3% methane in cO., mixture 
EOgo ss Cuce lia for other systems 


The temperature of the gas in the calorimeter is measured with 
an error of less than THOt01° Clitites AT = 0.01°C. The test value ¢ are 


related to the pressure a = P, -— (AP/2) and the inaccuracy of 


L 


determining the pressure is taken as +1.0 psi (the resolution of the 


i 1 (oe: a 
Heise gauge). The term , [Gp rx AP] = €) will be, 
+0.0003 a | < 0.002% for nitrogen 
and +0.0005 oees < 0.15% for 42.32 CH, in CO, mixture. The 


term | ie 


) AT ] = €,, will be less 0.1% for methane-carbon dioxide 
eee m 2 


mixtures. 
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APPENDIX B 


ENTHALPY FORMULAE AND LOOP CHECKS 
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I. FORMULAE USED IN CALCULATION OF PROPERTIES FROM EQUATIONS OF STATE 


lL. Benedict-Webb-Rubin Equation 


a) Enthalpy Departure 


4C4 Bane Gan0n 
(Hp — Body = (BoRT - 2A - gil el RT Ge ou) Peete 
2 -p? “yo? 2 
e ee oe 
pO PRS By yet 
T ye 


b) Isothermal Throttling Coefftetent 


2 
4C ry (2 
0 4 ce 3 250 
Re (ByRT-2A)- ay) + (2bRT-3a)p + bap” + —py— [5et5ye -2y PJ 
7 2 2.6 
RT + 20 (B.RT-A,- ls + 30 7 (bRT-a) + Rion fg ete [30 2 (1-+yp ) - 2y°e ] 
0 0 2 72 
2. Soave-Redlich-Kwong Equation 
a) kEnthalpy Departure 
T & - a Vib 
(ire AE Ne ; 1n 


b) Isothermal Throttling Coeffictent 
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II. THERMODYNAMIC CONSISTENCY CHECKS 
The isothermal data of two mixtures of methane and carbon 
dioxide of this work were used in conjunction with isobaric data of 
Bishnoi to test the consistency of the experimental data. The location of 
the loops is shown in Figure 34 and 35. The enthalpy change along 
different paths are given in Table XXVII and XXVIII. The maximum deviation 


was 0.37 percent. 
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A-28 
TABLE XXVII 
CONSISTENCY CHECKS ON THE MIXTURE OF 14.5 
MOLE PERCENT METHANE IN CARBON DIOXIDE 
Loops AH J Sain % Deviation 
1- 2-9 +2464.3 
1-10- 9 -2446.3 +0 .37 
2- 3- 8 +2440.7 
2- 9- 8 -2451.3 -0.21 
3- 4- 7 +3252.1 
3- 8- 7 -3243.5 +0.13 
4- 5- 6 +3938.5 
4- 7]- 6 —3924.7 +0 .18 
10- 9-12 +2264 .6 
10-11-12 -2261.9 +0 .06 
9- 8-13 +2327.4 
9-12-13 -2336.6 -0.20 
8- 7-14 +3867 .3 


8-13-14 =385940 +0 .36 


A=29 


TABLE XXVIII 


CONSISTENCY CHECKS ON THE MIXTURE OF 42.3 


MOLE PERCENT METHANE IN CARBON DIOXIDE 


AH J ole % Deviation 

+1943 .9 

-1935.5 he 
+2160.3 

-2150.4 +0 .23 
+2389.8 

-2404.4 -0.30 
+2617.1 

-2611.3 +0.11 
+1654.6 

-1646.4 +0.25 
+1902.1 

-1905.4 -0.09 
+2214.7 

-2214.4 +0.007 
+2450.1 

-2453.0 -0.005 
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TABLE xoxT 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR NITROGEN 


TEMPERATURE (K) = 273015 
PI AP POWER FLOW - AH -(AH/ WP) 
PSIA PSI J/MIN MOL /MIN J/MOL J/MOL-BAR 
42665 176.6 190.115 2 «1936 86.666 aap Rls 
607.5 201.4 282.256 229607 95 2334 6.865 
613.4 20 3.4 290.666 2.9842 97.400 6.945 
B13.4 205.0 339.198 3.6109 93 .936 6 2646 
B13.5 203.5 331.240 3.5841 92.419 6.586 
1112.0 PAGS: 401.705 43842 91.625 6.253 
1117.0 207.5 393.090 463947 89.445 6.252 
1415.0 211.5 432.904 5.0767 85.272 5.847 
1416.0 207.3 4126422 4.9977 82.520 5.773 
1416.0 20 6 0 410 «424 5.0208 81.743 5.755 
1713.0 178.2 342.220 5.2161 65.607 5.339 
1713.0 AS 338.782 5.2024 65.119 5.321 
1713.0 176.2 333.636 5 1748 64.473 5.307 
20 20 .0 147.6 249.642 5.1748 48.241 4.740 
20 20 «0 146.5 248.226 5 1470 48.227 4.774 
20 20 20 145.0 242.927 5.1190 47455 4.746 
2317.0 132.5 198.180 5 1748 38.297 4,192 
431.5 173.0 189.505 2 61964 86.276 7.233 
431.5 172.9 187.716 262049 85.133 7.141 
431.5 leieed 188.537 2 1964 85.836 7.208 
431.5 1i2e6 191.484 222106 B6 2620 7.278 
431.5 172.6 186.866 2.2021 84.856 7.130 
618.5 188.0 259.678 229374 88 403 6.820 
618.5 187.0 255.522 2 «9327 87.2127 6.757 
618.5 184.9 DoPeclt 2.9163 86 483 6.783 
813.5 186.5 293.902 3.4953 84.085 6.539 
813.5 184.7 289.751 304748 83.385 6.547 
1107.0 185.0 328.327 4731 78.675 6.168 
1107.0 183.8 326.156 4.1731 78.155 6.167 
1414.0 186.3 355.320 4 8055 73.938 5.756 
1714.0 186.2 557.066 5.3453 104.214 Sey 
1714.0 185.8 356.000 5 03453 66.599 5.198 
1714.0 185.7 357.251 5.3453 66 6833 5.219 
Neh Ate 185.5 352.784 5 63320 66.163 5.173 
20 14.0 164.9 293.360 5 64640 53.689 40722 
2014.0 164.0 289.930 5 4588 eis ili 4.697 
20 14.0 162.6 287.597 54301 52.963 4.724 
230 4.0 143.7 225.535 5 4562 41.335 4.172 
2304.0 143.0 D220 54379 40.846 4.142 
230 4.0 142.3 220.159 5 04248 40.583 4.136 
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TABLE XXXI (continued) 


ING hs. 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR NITROGEN 
TEMPERATURE (K) 


PI 


PSIA 


415.3 
425.3 
463.5 
463.25 
815.3 
615.3 
616.0 
619.5 
619.5 
GUT ew 
815.0 
817.0 
1115.0 
RL ae: 
1412.0 
1415.3 
LAL. 
VALS 
20 18.4 
20 21.3 
235 565 
241364 


AP 
PST 


178.0 
182.0 
1 eKehes tL 
197.0 
2Uo6 1 
208.0 
249.20 
211.5 
20 7 20 
206.9 
ELSI 8) 
2 20 t 
204.8 
NESS 
2256) 
206.4 
198.7 
183.5 
149.2 
152.3 
UZ Ziel 
107.7 


POWER 
J/MIN 


s9s001 
149.643 
178.370 
175.081 
213.546 
225.506 
29220 
230.857 
223.696 
222.568 
301. 703 
305.814 
eh SUS 
304.097 
348.631 
308.324 
292.300 
264.414 
192.258 
197.369 
135.247 
108.010 


BOS elo 


FLOW 
MOL/MIN 


200399 
2-2-0838 
2 62331 
2202998 
3 5260 
228728 
Sie 245 
229280 
2 8093 
228093 
36452 
32-6784 
4.1685 
422306 
4.9740 
4.7696 
SS Oy 
520426 
4.9691 
4.9897 
4.8283 
425983 


AH 


J/MOL 


68.576 
1.810 
79.874 
76.127 
TPS Les 
78 6496 
93.493 
78.842 
Use Ae) 
19.2224 
82-1766 
Beale 
69.034 
71.879 
70.089 
64 642 
Biles WS 
52436 
38.690 
39% DOs 
28.011 
232489 


—(AH/ AP) 
J/MOL-BAR 


eyelets 
52122 
5.845 
5.604 
Det29 
Drala 
52445 
5.406 
be DUD 
Die DO 
Se, a0) 
Sresloit 
4.887 
4.2906 
4.548 
4.2542 
4.172 
4-144 
3.759 
3.766 
3.309 
3.161 
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TABLE XXxXI (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR NITROGEN 


TEMPERATURE (K) = 343245 
Pil AP POWER FLOW AH = (AAA) 
PSIA PSI J/MIN MOL/MIN J/MOL J/MOL-BAR 
433.6 T9562 116.445 220103 Biles oe a 4.302 
437.7 194.2 115.902 2 00282 572143 4.266 
437.7 193.7 115.608 200282 56.998 4.266 
437.8 197.2 149,322 2.0521 58.145 4.276 
618.6 205.4 1153 2 4D4 226657 572564 4.064 
618.8 213.5 163.650 Ziel 259 60.092 4.082 
619.7 213.0 161.977 227095 59.780 4.070 
619 i 2265 160.043 2/095 59.066 4.031 
619.7 2l2et 160.764 267095 59.032 4.045 
619.7 212.7 161.977 227095 59.780 4.076 
814.8 2a 202 184.680 362676 56.517 3.862 
821.7 215.2 189.234 3.3061 57.2236 3-856 
818.6 225.26 201.634 323716 59.803 32843 
821.7 214.7 L8i7 .693 3.2954 56.955 32846 
1117.6 219.0 215.542 4.0098 ake Peel SYS) 32558 
1419.6 190.3 185.521 4.2858 43.287 Be299 
1723.5 24.6 of2 230.575 5.0837 45.355 32042 
20/1545 170 .0 157.410 4.9053 322089 Ze St 
2330.5 154.0 133.093 520421 26.2729 22497 
2410.5 116.8 85.158 4.4146 19.289 2-394 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR NITROGEN 
TEMPERATURE (K) = 374.215 


Bit AP POWER FLOW AH -(AH/ AP) 
PSTIA Pool J/MIN MOL/MIN J/MOL J/MOL-BAR 

e 06.3 97.202 19710 49.315 32467 
rents atau Ole oe 1.9557 46.2904 3-401 
438.8 199.5 90 . 146 We 955 46.398 Sep ete} 
438.8 199.5 89.227 eID t 45.622 cy ohio} 
623.8 210.5 120.692 205715 46.933 Benes 
819.8 210 0 136.690 3.0971 44.133 3.047 
1117.8 198.8 140.176 346508 Shabrjedel® 2-800 
1427.8 20461 153636! 4 22426 36.2149 2-568 
1717.8 210.8 163.650 4.1643 34.349 2-363 
1717.8 199.3 161.677 4.7519 34.023 2.475 
20 13-8 194.5 141.971 5.0038 28.312 Zetls 
162.1 194.853 4 8254 21./29 1.944 
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A-37 
TABLE XXXII 
TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
THE 14.5 MOLE PERCENT OF CH4 IN CO2 
TEMPERATURE (K) = 363.15 

Pl AP POWER FLOW AH —( AH/ AP) 
PSIA PST J/MIN MOL/MIN J/MOL J/MOL-BAR 
463.6 194.5 599 429% 1.7574 340.996 254427 
515.6 20 362 670.778 1 8832 356.175 254416 
611.6 197.0 768.228 201577 356.024 26-211 
611.6 19652 Oe 6097 221497 354.778 26.219 
613.6 181.7 664.204 220273 327.2616 262144 
809.6 175.6 825.870 2 4907 331.578 27.386 
813.6 182.5 845.143 204602 343.519 27.2300 
1912.0 174.2 994.177 28914 343.837 28.619 
1012.0 172.5 946.085 221823 340.035 28.590 
1212.0 162.6 1957.043 321641 334.067 29.798 
1216.0 15055 979.680 320210 324.286 29-862 
1412.0 158.2 1165.951 3 4688 336.118 30.815 
1406.0 138.7 928.892 321363 296.173 BO.959 
1406.0 PS0r a 847.899 3.0333 279-527 31.090 
1612.0 145.22 1158.108 306477 317.483 31.712 
1612.0 141.0 1199.409 505910 308.880 Shet he 
1612.0 he) 567.449 Zot 53 2032146 32.2201 
1818.0 140 .0 1211.007 329045 310.154 32.131 
1994.0 118.0 1001.584 328234 261.961 32.198 
2) 20 .0 M2Ss2 1098.551 39755 276.330 Bh.996 
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A-38 
TABLE XXXII (continued) 
TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
THE 14.5 MOLE PERCENT OF CH4 IN CO2 
TEMPERATURE (K) = SSeS 
PJ AP POWER FLOW AH -(AH/ AP) 
PSIA PSI J/MIN MOL /MIN J/MOL J/MOL-BAR 
415.6 169.8 692.946 1.6171 428.501 36-601 
605.6 184.5 1134.688 222090 513.666 40.380 
605.6 184.0 1130-029 222066 512.100 40.366 
813.6 163.5 1368.037 226132 523.2509 46 6439 
1910.0 143.4 1578.056 229236 539.746 54.591 
10 10.0 142.5 1567.262 209155 537.552 54.712 
1214.0 134.0 2047.897 3.3280 615.335 66.602 
1214.0 133.2 20 31.559 323210 611-727 66.584 
1310.0 128.8 2295-553 325142 653.212 3.556 
1412.0 12763 2671.699 327815 706.508 80.495 
1412.0 126.5 2647.914 32/7736 701.683 80.451 
1415.0 83.0 1450 .356 32-0791 471.026 82.309 
1524.0 106.1 2314. 156 307154 623.000 85.163 
1524.0 104.9 2279.477 323-7003 616.013 85 oll 
1611.0 115.6 2776. 365 4.1741 665-4140 832452 
1701.0 109.2 2422.618 4.23113 561.915 14.2632 
170 1.0 108.2 2378.256 4.2854 554.960 74.390 
1814.0 96.0 1699.500 422540 399 .502 60.357 
1915.0 102.9 1577.787 4.5818 344.354 48.536 
1915.0 102.5 1570-591 4.5735 343 .410 48.592 
1915.0 101.2 1536.208 4.5553 S33 (253 48.331 
20 13.0 86.0 981.341 4.2963 228 0412 36 .5:21 
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TABLE XXXII (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 


a | 
PSIA 


513.6 

71326 

915.6 

915.6 
1112.0 
1112.0 
1314.0 
1314.0 
1314.0 
1512.0 
1514.0 
1514.0 
1762.0 
1762.0 
1762.0 
2012.0 
2012.0 


POWER 
J/MIN 


931.854 
981.050 
10 73.432 
10 65.170 
1376. 845 
1372.038 
1639.549 
1635.234 
1626.892 
1660 .767 
1798.145 
1769.884 
1770 .634 
180 2.599 
1730. 858 
1329.178 
1290.213 


Serio 


FLOW 


MOL /MIN 


1.9674 
23166 
206157 
226118 
3 20448 
320414 
34118 
344099 
34010 
325912 
327095 
326826 
3 29594 
329817 
3.9369 
329817 
329427 


THE 14.5 MOLE PERCENT OF CH4 IN C02 
TEMPERATURE (K) 


AH 


J/MOL 


4732635 
423.471 
410.375 
407.829 
452.2192 
451.107 
480.538 
479.550 
478.345 
462 6452 
484.129 
480.600 
447,197 
452.714 
439.649 
333% s8 V6 
327-236 


A-39 


—(AH/ AP) 


J/MOL-BAR 


Bho BIT, 
34.602 
37.814 
Dre V6 
41.196 
41.175 
45.2140 
45.164 
45.197 
49.137 
48.822 
48.950 
51.071 
SV 291 
51.012 
47.097 
47.131 


TABLE xXxxII (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
THE 14.5 MOLE PERCENT OF CH4 IN Ch2 
TEMPERATURE (K) 


Pil 
PSIA 


418.5 
418.5 
419.5 
420.5 
Bi 65 
ble 5 
555.65 
609.5 
613.5 
708.5 
807.5 
957.6 
O50 «9 
959.5 
959.5 
1009.6 
1914.0 
1054.0 
1052.0 
1054.0 
1108.0 
1108.0 
1108.0 
1156.0 
120 4.0 
1204.0 
1234.0 
1262.0 
1262.0 
1298.0 
1318.0 
1318.0 
1318.0 
1318.0 
1362.0 
1614.0 
1614.0 
1614.0 
1464.0 
1813.0 
1813.0 
2008.0 
2009.0 


AP 


POWER 
J/MIN 


830.148 

825.384 
823.012 
840.886 
19 36.280 
10 27.045 
19 28.620 
1142.614 
1158.057 
1470 .338 
180 7.548 
1827.272 
2817.071 
2779.054 
2753.518 
22432140 
2359.6/0 
23022114 
226% e335 
2250 .984 
3348. 740 
3875.93) 
3362.466 
4276.152 
4685.280 
4551.967 
5017.247 
40 43.094 
409 9.822 
3654.539 
30 54.063 
30 13.848 
3000.119 
2933 3299 
2367.385 
951.568 
939.028 
920.727 
1298.474 
566.047 
566.047 

311.39¢ 
293.467 


293%,15 


FLOW 
MOL /MIN 


1.2 
1.7067 
1.6562 
1.6770 
1.9754 
1.9631 
2 -0659 
Zeb 992 
201467 
264938 
2 06698 
1.7617 
3.0885 
321005 
321081 
1.9000 
2 69580 
220906 
204770 
224568 
28421 
208942 
2etih 
320768 
3.2941 
Se296F 
3 8250 
328155 
3 8311 
421202 
4.1635 
4.1525 
4.1399 
4.1051 
4.3587 
5.0363 
5 .0300 
5.0000 
4 25320 
5.0858 
5.0889 
4.217829 
4.6643 


AH 
J/MOL 


485.108 
483.607 
496.904 
501.398 
524.579 
5123 01 ST 
497.883 
519.543 
539.437 
589.575 
677.2030 
1037.166 
912.089 
896 .305 
885.899 
1180.577 
OK ZL 
1101.130 
OVSis2 42 
916.193 
1178.242 
1166 6423 
1168.420 
1389.764 
1422.285 
1380.748 
1311.665 
1059.625 
1046.641 
886.974 
1334520 
725.776 
724.667 
714.607 
543.137 
188 .941 
186.682 
184.144 
286.509 
Die 293 
Dus 231 

65.105 

62.2916 


A-40 


—( AH/ AP) 


J/MOL-BAR 


44.871 
44.904 
45.327 
45-224 
48.307 
48.608 
50.675 
53.823 
Doerr 
60.007 
68.908 
BOR eto 
Vile 930 
92.2001 
91.843 
103.774 
1735209 
124.770 
130.668 
130% 022 
166.884 
166.675 
167.2455 
2154587 
236.837 
23s 216 
2eteet0 
190.914 
189.753 
152.062 
1286959 
128.372 
128.646 
IZSetot 
935115 
28 «634 
28-441 
28 6412 
50.065 
17.451 
17.440 
11.952 
let 


E57 


FEN. TSI 
1H 28 LS 
VERLIES 
A¥S 2 TES 
OS .LS8 
afe.0eT 
fay eet 
sab, Sel 
ade’, 3S] 
SVE. BSI 
OP ae BSE 
fav. ht 
At; s€0 


yaalees 


O.A@f af 
0.°f af 
GO, fal 
0, dei 
O,61 aL 
O.€ fa! 
0.8005 
0,2 005 


TABLE XXXII (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
THE 14.5 MOLE PERCENT OF CH4 IN CO2 
TEMPERATURE (K) 


PY 
PSIA 


419.5 
419.5 
483.5 
5D.) 
628.5 
629.5 
O95 si? 
745.25 
73 9D 
T7B65 
810.5 
1246.5 
1246.5 
1265.5 
1265.5 
1275.5 
17293615 
1293.5 
1308.5 
1855.5 
1414.5 
1473.5 
1473.5 
1529.5 
1595.5 
1595.5 
1655.5 
1705.5 
IGS 65 
La15 65 
1815.2 
1909.5 
2009.25 


AP 


PSI 


POWER 
J/MIN 
889.046 


887.201 
10 36.845 


191.556 


1366.217 


1367.535 


1529.424 


1772.513 


1955.87 7 


1944.557 


2129.922 


1219.439 


993.305 
919.281 
852.661 
839.516 
821.191 
896.192 
739.603 
591.679 
506.595 
421.487 
410.706 
380.994 
320% 825 
317.986 
266.851 
254.394 
2500's 192 
223.091 
218.274 
166.816 
126.672 


283.15 


FLOW 
MOL/MIN 


1.6090 
1.6107 
1.7854 
1.9573 
221245 
201245 
22476 
2 o3803 
224638 
2 04596 
Ze D219 
5 4257 
bot 261 
5 1564 
5 0463 
5 .0833 
521534 
5.1504 
5.1017 
5 20245 
50864 
5 0419 
5.0026 
5 00864 
5.0431 
5.0214 
4.29362 
5.0182 
429806 
5 0956 
Bie Oi ie 
4.8295 
426552 


AH 
J/MOL 


5520026 
550.815 
580.727 
608.755 
643.058 
643.679 
680 -442 
744.2628 
793 815 
790.571 
844.54) 
224.750 
193 oft) 
178.276 
168.967 
165.149 
159.347 
156.528 
144.969 
127 .a58 
99.207 
83.596 
82.098 
74.903 
63.615 
632326 
54.060 
50.693 
Die 255 
43./80 
42.991 
34.541 
27.2210 


A-41 


-(AH/ AP) 


J/MOL-BAR 


54.963 
54.868 
58.654 
63.749 
2.896 
71.265 
80.235 
89.182 
96-104 
96.113 
108 .687 
32.435 
31.226 
28.666 
28.331 
27.343 
25-822 
25.508 
24.112 
202454 
17.196 
14.822 
14.791 
13.248 
14°. 562 
14.596 
10.316 

9.462 

9.536 

8.068 

TLCREM Ae: 

72146 

6.118 


158 


TABLE XXXII (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
ThE 414.5 (MOLE PERCENT (GF (Ch4 TIN (Co? 
TEMPERATURE (K) 


PI 
PSIA 


413.5 
413.5 
478.5 
DZ Die 
Dililie'D 
585.8 
1051.0 
1051.0 
1100.0 
1198.0 
1228.0 
1256.0 
1311.0 
1354.0 
1398.0 
1398.0 
1442.0 
1526.0 
1594.0 
1657.0 
1657.0 
1711.0 
1711.0 
1711.0 
i ooe0 
IMT r aw 
1816.0 


/AP 
Pin 


138.0 
S75 
141.6 
3260 
12D ist 
122.6 
89.5 
82.7 
80.5 
80 .9 
19.4 
82.5 
78.8 
80.7 
Wal 3D 
Til 09 
1305 
78.8 
8125 
78.2 
7725 
76.5 
eo.) 
74.6 
Uaioe 
1320 
66.2 


POWER 
J/MIN 
DD 28 


941.973 
1214.074 


1392.969 


1424.137 


1464. 387 


1278.791 
991. BOD 
343.843 
282.580 
256.801 
266.048 
216. 733 
ZUG. 5D 
185.859 
188.255 
161.808 
156.197 
156.197 
133.611 
132.084 
118.586 
117.876 
112.882 
103.790 
191.763 

78.350 


273.15 


FLOW 
MOL/MIN 


1.6319 
1.6254 
1.8465 
1.9452 
220389 
2.0733 
522261 
5 20250 
5.0062 
5.0623 
5.20125 
5.1299 
4.9365 
520932 
4.9999 
5.0125 
4.8786 
4.9873 
D tli 
5 20809 
520561 
Src 0204 
5.0094 
4.9684 
429410 
4.9173 
4.26995 


AH 
J/MOL 


583.803 
51 9%,D09 
657 495 
669.829 
698 477 
706.282 
244 .688 
198.575 
68 .682 
55.819 
Die cot 
51.862 
43.903 
43.048 
37.2172 
BY sO 20 
33.166 
B41 .308 
30.170 
262296 
262123 
232584 
23505") 
A Pal 
21.005 
20.694 
16.671 


A-42 


mo AtAl/ APA) 


J/MOL~BAR 


61.357 
61.127 
67.345 
Chek poe § 
80.850 
83.554 
39.6512 
34.825 
12.374 
10.007 
Mie S58 
SP te yf 
8.080 
7736 
6.956 
6.992 
62544 
5.164 
5-369 
4.877 
4.888 
4.471 
4.496 
4417 
4139 
4.111 
32652 
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FABLE XOOCriT 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
THE 42.3 MOLE PERCENT OF CH4 IN CO2 
TEMPERATURE (K) 


Pe 
PSIA 


431.4 
431.4 
429.4 
619.4 
815.4 
815.4 
1917.0 
1017.0 
1217.0 
1217.0 
1217.0 
1415.0 
1613.0 
1613.0 
1811.0 
1811.0 
2017.0 
2017.0 
2017.0 


AP 
PSI 


183.8 
186.2 
183.0 
195.0 
176.8 
106.2 
167.5 
166.8 
168.2 
167.5 
167.0 
161.4 
164.5 
164.0 
148.8 
148.0 
126.9 
121.8 
119.6 


POWER 
J/MIN 


434.626 
431.792 
429.602 
616.525 
653.316 
649.560 
199.719 
705.380 
812.023 
806.890 
89 1.588 
842.537 
939.909 
936.676 
859.701 
857.203 
1790 7.543 
668.143 
649.549 


263615 


FLOW 
MOL /MIN 


1.8047 
1.7982 
1.8017 
223833 
21407 
20 1356 
3.0831 
3-2-0761 
34885 
324760 
304677 
327612 
4.1129 
4.1094 
4.2041 
4.1971 
4.1517 
4.0738 
4.0197 


AH 
J/MOL 


240.829 
240.114 
238-442 
258 .683 
2382370 
2372442 
250.492 
229 2302 
232.1768 
232.128 
284.257 
224.007 
2282522 
Lee I52 
204.508 
204 2233 
170.418 
164.009 
161.589 


-(AH/ ®) 


J/MOL-BAR 


19.004 
19.009 
18.897 
19.240 
19.554 
19.544 
19.952 
19.938 
20.071 
20.099 
20.075 
20.129 
20.148 
20.157 
195932 
20.014 
19.477 
Ware 30 
19.595 


uw 


— aie yp * 


TABLE XXXIII (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
THE 42.3 MOLE PERCENT OF CH4 IN CO2 
TEMPERATURE (K) 


PI 


PSIA 


426.5 

619.5 

820.5 

820.5 
1015.0 
1215.0 
1215.0 
1215.0 
1215.0 
1415.0 
1415.0 
1613.0 
1613.0 
1813.0 
1813.0 
2915.0 
2915.0 


AP 


PS 


LEB e2 
190.5 
184.2 
183.6 
170 .8 
169.7 
169.2 
168.3 
UGifkeré 
167.0 
166.0 
150.2 
149.5 
152.3 
150.5 
128.0 
117.4 


POWER 
J/MIN 


59 6.400 
781.843 
936.376 
931.765 
109 4.805 
1154.349 
1150 .465 
1140 .957 
LBZ Fish DD 
1276.662 
1265.803 
120 2.488 
1190.880 
1309.613 
1288. 312 
1045. 762 
De 316 


5S Ble kD 


FLOW 


MOL /MIN 


1.8583 
205147 
269990 
229956 
3 23490 
Ba lt(l 
3.7661 
367622 
321544 
421448 
4 1360 
422988 
4 «2886 
4.26685 
4 26405 
4.5934 
44027 


AH 
J/MOL 


2724492 
310.903 
BUD eZ26 
SI S035 
300.031 
305.565 
305.476 
303 .264 
300.854 
308 .008 
306.040 
279.2124 
277.683 
28.0 6° 187 
277.620 
227.2663 
207.443 


=CAHT 2?) 


J/MOL-BAR 


22.818 
23.670 
24.584 
24.570 
25.477 
26 #f 5 
26.2185 
262134 
ZOreST 
26.750 
2Oera9 
27-011 
266939 
26.714 
26.1754 
25.796 
252627 


SEStEeis & 


: 


Sbaneets 


A-45 


TABLE XXXIII (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
THE 42.3 MOLE PERCENT OF CH4 IN CO2 


TEMPERATURE (K) = 21 3b 
PI AP POWER FLOW AH =(AH/ Pe) 
PSIA Sy] J/MIN MOL /MIN J/MOL J/MOL-BAR 
427.4 166.5 583.498 e290 302.479 264348 
427.4 166.0 BM feds hohe, 1.9224 301.248 26 320 
427.4 165.8 5116s 903 1.9224 300.605 26.296 
4274 Voor 2 574.705 1.9158 AS}S/ BST 26.336 
61524 183.0 99 6.149 2 5860 BOoO.395 272770 
615.4 Tore6 902.768 225807 349.807 27.184 
615.4 182.5 899.968 2257154 3492442 Zi ongia 
813.4 N60 979.234 229702 329.680 29 De?o 
813.4 161.4 978, 313 20 F102 329.370 29,298 
1017.0 148.0 10 70 .568 323488 BLS Ou 31 ei. 
1217.0 145.9 1253.064 32/766 S33 eto? 326905 
1417.0 144.3 1431. 353 4.1623 343 884 34.564 
1617.0 134.0 1435.355 4 4429 323.060 34.967 
1617.0 pRB 1402.092 424023 318.488 34.941 
1715.0 126.8 1367.427 4.5110 303.128 34.672 
1715.0 126.0 1356. 21 2 4.4917 301.949 34.757 
1715.0 124.5 1330.150 4 04658 291.822 34.698 
1815.0 127.4 1403.239 4.7025 2 IB e399 Saco 
1815.0 126.8 1392-6201 4.6902 296.828 eloysleys 
1815.0 126.4 13836312 4.6871 2 Sete O 33.864 
1815.0 124.3 IBS Zee 4 26500 290.904 335943 
20 21.0 103.7 10 28.901 425874 2242285 31.369 
2021.0 LO? 1 1010.614 4.5621 Ze leo 1. Bate oct 


20 21.0 102.0 100 3.029 425440 Z20eloS only Sie) 7 


1% VHA = 7 IMC 


abs. S08 
eae, 108 


OPE. HPS 
BSB. HFS 
asf, ces 
ae. GOS 
2894668 
£S@s £88 
EET.O8S$ 


TABLE XXXIII (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
THE 42.3 MOLE PERCENT OF CH4 IN CO2 
TEMPERATURE (K) 


PI 


PSIA 


424.6 
424.6 
615.6 
815.6 
815.6 
10 16.0 
1916.0 
10 16.0 
1118.0 
1118.0 
1118.0 
1214.0 
1214.0 
1214.0 
$322.0 
1322.0 
1322.0 
1416.0 
1418.0 
1418.0 
1514.0 
1514.0 
1612.0 
1612.0 
1714.0 
1812.0 
1812.0 
1812.0 
1914.0 
1914.0 
20 16.0 


AP 
re 


147.5 
147.3 
15% 62 
157.0 
PS fen 
148.5 
147.7 
147.5 
134.4 
143.3 
18265 
126.5 
124.5 
125.4 
114.8 
Lise 
123.0 
110.0 
109.3 
108.7 
10.135 
101.2 
111.0 
110.3 
99.0 
92.6 
926 1 
90 .4 
85.8 
80 .0 
7626 


POWER 
J/MIN 


791.026 
789.374 
130 7.205 
1960 .814 
1968.95 
2782.203 
20 U2. L5G 
2760 .099 
30 26.130 
2993.953 
2971.018 
3331. 837 
3266.655 
3290 .013 
3158.507 
ST 286200 
SiS 921 
2195.069 
2777.665 
2757.898 
2117. 643 
2106.894 
1998. /08 
1972.888 
bS2peor4 
973.487 
sheyeleps hots} 
S10) ry tlleye: 
696.325 
617.927 
477.233 


273.15 


FLOW 
MOL/MIN 


2.0285 
2.0285 
2e74e1 
Bala 
3.4177 
4.0433 
4 20266 
4.0229 
4 1836 
4.1766 
421766 
4.3918 
4.3751 
403751 
4 65260 
4.5356 
405356 
407454 
407454 
4.7485 
4 08746 
4.8717 
5 3951 
5.3714 
5 63501 
5 03736 
5 63603 
5.2828 
5 2867 
5.0575 
5.0143 


AH 
J/MOL 


339.4936 
389.124 
475.668 
573.413 
576.100 
688 .090 
688.457 
686 .083 
£23..885 
716.826 
TUY.333 
758.640 
746.632 
TORS ILT 
697.853 
689./00 
686.098 
588.997 
585.330 
580.791 
434.415 
432.473 
370.462 
3672289 
2472836 
NG1eiko9 
179.092 
176.068 
13r.710 
phedrdes shies) 
95.173 


A-46 


—=AAta/ oy 


J/MOL~-BAR 


38.342 
38.314 
43.886 
53-000 
53.05% 
67.204 
67.604 
67.463 
78.056 
172994 
17.864 
86.981 
S6e979 
86.913 
88.166 
88.134 
88.062 
77.660 
772671 
TT 2494 
62.075 
61.981 
48.406 
48 2296 
36.308 
28.374 
28.203 
28 2248 
22-26% 
2oei OU 
18.020 


163 


0% \Hn 
AAG SHANE 


Soe. RE 
xb? ak 
ye 


f oes 
o7s.0e39 
Met a 

we. To 


ries 
dhe ST 
piss 48 
CTS .38 
eyn,43 
OAtE BS 
aé& [4 8A 
SAU eas 
Osa. TV 
(Vv Aaeyt 
ea, FT 
270.509 
rae ste 
d0e 68S 
4PS, 8h 
RIE «4% 
ere. 8s 
£08.85 
Pel, Rs 
OaS 25 
QefeSs 
080,85 


ae 


AEB, FHS 
PAT.IRS 
$e0.eTi 
Gao, ay i 
OFV.hE! 
HVIeSst 
Evi .e? 


A-47 


TABLE XXXIV 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
EQUIMOLAL MIXTURE OF CO2-CH4-N2 


TEMPERATURE (K) = 283.15 

PI AP POWER FLOW AH SAngnr 
PSIA PSI J/MIN MOL/MIN J/MOL J/MOL-BAR 
416.4 199%,5 441.619 i, o182 273/20 19.2899 
419.4 200 20 450.868 1.6249 2772459 206121 
52264 209.0 574.881 12-9451 295%909 20.510 
607.4 ZAZs5 667, 132 201783 306.529 206921 
121.4 192.0 661.506 223471 281.835 2167290 
806.4 13269 672. 163 2 «4837 270.866 21.526 
1009.4 172.7 740.144 228157 262.857 22.075 
1009.4 171.5 %314693 2 68108 260.311 222014 
1208.4 170 .4 826.584 321564 261.869 22.299 
1410.4 163.5 858.777 34302 2506355 22.208 
1410.4 163.2 855.692 324235 249.941 22.212 
1410.4 162.5 849.331 34182 248.470 22-177 
1558.4 163.2 897.053 326563 245 2341 214803 
1558.4 162.2 890.252 3.6451 244.2229 21.838 
1715.4 165.5 943.1792 309134 241.2163 21.134 
1725.4 165.0 938.749 329010 240.016 rag WO he De | 
1858.4 164.0 938.985 421022 228.897 20-243 
1858.4 163.2 935.717 4 0934 228.2589 20.315 
20 29 64 146.1 7772128 4.0857 190.204 18.882 


TABLE XXXIV (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 


PI 
PSIA 


415.26 
5636 
759.6 
715926 
908.6 
908.6 
1917.6 
1017.6 
1059.6 
1059.6 
1205.6 
1205.6 
1309.6 
8309.6 
1363.6 
1363.6 
1506.6 
150 6-6 
1506.6 
1658.3 
1658.3 
177926 
1909.6 
201326 
2013-6 


TEMPERATURE (K) 


AP 
PSI 


194.2 
216.0 
208.0 
20 6.0 
193.7 
191.4 
176.6 
175.6 
194.8 
192.0 
19562 
190.8 
170 .0 
169.2 
183.0 
181.0 
15265 
149.6 
146.3 
149.6 
148.8 
ier) 
Liz en 
136.3 
135.1 


POWER 
J/MIN 


541.518 
7995529 
950.825 
945.096 
BB es OS 
975.728 
10 56.383 
1051.193 
12236240 
119 7.195 
1222.936 
1194.714 
1244. 736 
1229.464 
1200 .628 
1183. 325 
10 36.447 
119 8.530 
1016.5/2 
1128.871 
1118.598 
109 9.040 
918.665 
896.962 
882.717 


EQUIMOLAL MIXTURE OF CO2-CH4-N2 
263.15 


FLOW 
MOL/MIN 


1.6735 
2ell21 
224959 
2 04997 
226657 
26761 
320580 
3 20550 
22-9281 
22-9312 
321178 
321223 
326045 
Seog 
322596 
32610 
304354 
321498 
Ze 2bC2 
4.0108 
309995 
4.0602 
4.1688 
4 23246 
43111 


AH 


J/MOL 


Sear a3 
378.546 
380.944 
378.083 
369.619 
364.595 
345.445 
344.088 
383.604 
Site 16 
B92 625) 
BO 26O34 
345.319 
3422167 
368.331 
362-866 
301.693 
295.621 
289.028 
281.452 
279.679 
248.518 
220.363 
207.407 
204.754 


A-48 


= (A th aPA 
J/MOL-BAR 


24.166 
25.2418 
26.563 
262619 
27-676 
272628 
28.370 
28.420 
28.561 
28.532 
29.098 
29-086 
29-461 
29-330 
29 6 2 
29.077 
28 .693 
28.660 
28.6293 
2724286 
272260 
252 16% 
24.2121 
22.070 
22 OL 


TABLE XXXIV (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 


PI 
PSTIA 


418.3 
513.3 
604.3 
TL2ed 
712.3 
863.3 
863.3 
100 4.3 
1004.3 
1109.3 
1109.3 
Paks ele e: 
1208.3 
1208.3 
MS (es nc 
1559.3 
1707.3 
1707.3 
1865.3 
1865.3 
20 16.3 


EQUIMOLAL MIXTURE OF CO2-CH4—-N2 


TEMPERATURE (K) 


AP 
PESst 


Tre 
192.4 
174.5 
179.4 
Ve D2 
152.5 
151.7 
156.0 
155.5 
158.1 
157.0 
156.3 
P5or #1 
156.7 
149.7 
142.3 
143.5 
142.5 
133.5 
see 
136.2 


POWER 
J/MIN 


665.239 
855.550 
881.086 
11034:229 
1103.721 
119 3.743 
1094.312 
L372.525 
1365.143 
1562.947 
155720912 
1537.906 
1689. 764 
1665.987 
1642.378 
1418.202 
1282.090 
1276.394 
Mira Vey: 
963.597 
850 .046 


243.215 


FLOW 
MOL/MIN 


1.7545 
20720 
202571 
205196 
225796 
rela 
20 1664 
321506 
301474 
3 24032 
324112 
3.3898 
325957 
3.5970 
Ze 299 


3.9931. 


42377 
4.2377 
422619 
4.2755 
404720 


AH 
J/MOL 


379.149 
412.899 
390.353 
427.667 
427.857 
398.009 
395.560 
435.632 
433.731 
459.249 
454.986 
453.682 
469.929 
463.154 
429.343 
355-160 
302.538 
301.194 
228.077 
225.373 
190.079 


(HARM) 
J/MOL-BAR 


29 644 
31.125 
322444 
34.575 
34.629 
372853 
37.818 
40.502 
40.454 
42.130 
42.032 
42.2099 
42.839 
42.868 
41.597 
362199 
30 O18 
30.655 
24.778 
24.2669 
20.241 


TABLE xXXxv 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
EQUIMOLAL MIXTURE OF CO2-CH4—-C2H6 
TEMPERATURE (K) 


PI 
PSIA 


270 4 
372.4 
426.4 
520 «4 
613.4 
113.4 
718.6 
903.6 
1107.6 
1214.6 
1214.6 
1236.6 
1253.26 
1407.6 
1413.6 
1417.6 
1413.6 
1458.6 
1611.6 
1618.6 
1613.6 
1657.6 
1814.6 
203166 


AP 


Pon 


146.3 


HO! Giet> 


211.2 
219.0 
215.6 


218.2 


195.2 
214.0 
215.8 
177.2 
175.7 
201.0 
191.0 
185.6 
191.0 
191.0 
172.7 
165.5 
164.6 
172.2 
150.8 
186.6 
129.4 
111.6 


POWER 


J/MIN 


"216.171 


417.219 
515.212 
6352649 
79 8.807 
8121/25 
699.629 
948.509 


1125.895 


905.836 
895.723 


1101.612 
10 26. 381 
10 71.682 
1123./14 
1122.845 


968.771 
256205 
968.889 


10 45.663 


854.476 


1182.802 


79 8.324 
577.2409 


© 363415 


FLOW 


MOL /MIN 


0.9227 
1.3179 
1.5023 
1.7711 
Le Oil Dit 
262088 
221125 
225710 
2.9639 
228686 
2.8591 
3.0901 
36030 
322446 
322943 
323011 
321396 
321432 
3 03434 
3643525 
32039 
326150 
3.2053 
321826 


AH 
J/MOL 


234.268 
316.092 
342.932 
358.885 
358.758 
367.935 
331.184 
368.917 
379.858 
315.767 
313.285 
356.493 
337 s695 
330.288 
341.105 
340.135 
308.558 
294 2347 
289.785 
304 -645 
266.698 
3272184 
2206979 
181.424 


A-50 


=My AP) 


J/MOL-BAR 


236224 
23.366 
23000 
23.168 
24.134 
242456 
24.607 
25.003 
25.530 
25.845 
25-861 
250i 23 
252640 
25.810 
252902 
25 2828 
ane 5's 
256195 
252534 
25.2659 
25-650 
25-430 
24.1768 
232578 


S¥e<Es 


TABLE XXXV (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
EQUIMOLAL MIXTURE OF CO2-CH4—-C2H6 
TEMPERATURE (K) 


ea | 
PSIA 


273-7 
Dee eth 
417.7 
515. 
D156 4 
1564 
715.7 
912.7 
912.7 
1113.7 
1319 6% 
1459.7 
1459.7 
1459.7 
1459.7 
1607.7 
1607.7 
1709.7 
L709 « | 
1709.7 
IBS. ft 
2019.7 
20:19.2 
20,19". ¢ 


AP 


PSI 


142.3 


164.6 
18% 6.5 
189.8 
189.6 
196.0 
196.0 
199.1 


197.3 


173.6 
164.6 
154.5 


eS es 


149.5 


142.8 


138.4 


139.5 


169.1 
164.2 
160.5 
142.8 
ji Me 
117.1 
117.1 


POWER 
J/MIN 


265.675 
375.110 
558.244 
677.809 
679.617 
942.049 
944.165 


120 7.128 


1192.003 


1193.59) 


1289.414 
1263.842 
1248.582 


1292.081 


1124.295 


Li3hb.o18 


1144.366 


1573.905 


1503.814 
1447.880 
1228.598 


999.537 


880.115 


880.115 


333.15 


FLOW 


MOL/MIN 


0.9754 
1.1702 
1.5061 
1.7642 
1.7616 
222560 
202566 
2/070 
226935 
2 «9265 
322200 
343635 
323525 
323131 
32355 
324114 
364285 
39348 
3.8719 
3.8371 
327628 
3 06942 
326551 
326526 


AH 


J/MOL 


272.357 
320,532 
370.651 
384.183 
385.785 
417.573 
418.396 
445.912 
442 2548 
407.816 
400.431 
375.771 
372 428 
362.824 
347.483 
331.686 
233.449 
3936994 
388 .389 
377.2336 
326.2506 
2464202 
240.788 
240.950 


yee th 


=(AH/ AP) 


J/MOL-BAR 


27 «59 
28.243 
28.671 
29.35% 
292518 
30.899 
30.960 


32.483 
32.932 
34.071 
35-284 
35.2215 
35.258 
35.199 
352292 
34.1759 
34.702 
34.307 
34.306 
34.098 
336 kG? 
29.881 
29.823 
29.843 


TABLE XXXV (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
EQUIMOLAL MIXTURE OF CO2-CH4-C2H6 
TEMPERATURE (K) 


PI 
ay | 


21305 
Bibeln) 
420.3 
Pp fad SO) 
Sled hg 
713.3 
918.3 
1113.3 
1215.3 
1313.3 
1313.3 
1413.3 
1413.3 
1413.3 
1414.3 
1414.3 
1511.3 
1511.3 
Lele 3 
1661.3 
1661.3 
179563 
1915.3 
1915.3 
20 28.3 


MP 


PSI 


129.4 
149.9 
171.6 
195.5 
195.1 
193.0 
194.0 
POeeL 
180.2 
168.0 
167.0 
157.2 
[Naito Pls) 
154.2 
146.5 
142.3 
143.0 
141.7 
141.4 
134.4 
133.1 
145.4 
139.3 
135.4 
126.3 


POWER 
J/MIN 


278.467 
390.678 
603.941 
873.902 
872.234 
1174.692 
1557.478 
199 4.919 
1945.115 
1910 .388 
1892.6 159 
1859.563 
1826.869 
189 3.783 
1672.413 
169 7.940 
1688 .840 
1668.011 
1652. 738 
1559.186 
1534.189 
1696, 760 
1478.298 
1412.976 
1160 .603 


Bisel 


FLOW 
MOL/MIN 


0.9723 
1.1670 
1.5211 
1.8788 
1.8788 
223573 
228598 
se 2926 
324353 
3 25354 
305253 
3 6338 
326153 
3.6017 
Ze 2072 
364613 
326558 
326411 
326356 
3.8199 
3.1954 
41944 
4.2869 
4.2277 
402225 


AH 
J/MOL 


286.385 
334.2763 
397.033 
465.115 
464.6227 
498.316 
544.593 
SB eis 2 
566 2204 
540.348 
536.729 
511.736 
50515305 
500.802 
476 .803 
464.544 
461.957 
458.096 
454.589 
408-166 
404.212 
404.526 
344.835 
334.212 
274.855 


A-52 


=—((Why AP? 
J/MOL~BAR 


32.099 
32.390 
33.557 
34.506 
34.510 
37.448 
40.714 
43.908 
45.572 
46.649 
46.614 
47.214 
476130 
47.104 
47.204 
47.348 
46 6854 
46.888 
46.628 
44.047 
44.046 
40.351 
35.903 
35.800 
31.563 


TABLE XXXV (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 


PI 


PSIA 


269.5 
215 oP 
365.5 
365.5 
419.5 
419.5 
613.5 
613.5 
815.5 
1917.5 
1017.5 
1112.5 
W125 
PUP2. 5 
121565 
oe ere) 
1415.5 
1613.5 
1613.5 
1613.5 
1817.5 
1817.5 
20 21.5 
2021.5 


TEMPERATURE (K) 


AP 
PST 


131.6 
152 12 


175.6 


175.0 


149 13 


178.6 


181.6 


181.0 


184.5 


186.1 


185.2 


178.9 


176.5 


174.5 


1572.7 


146.1 
127.6 
126.5 


126.3 


123.4 


119.4 


120.2 


117.8 


110.4 


POWER 
J/MIN 


344.346 
483.661 
676.531 
673.1792 
799.318 
789.560 
1222.887 
1214.028 
1801.682 
2589.961 
25762171 
2790 .391 
2768.562 
2741.524 
25872446 
2657.803 
2258.605 
1870.031 
1869.417 
1796.158 
1223.165 
12326599 
839.696 
761. 704 


FLOW 
MOL/MIN 


1.0122 
1.2128 
124485 
1.4437 
1.6193 
1.6164 
221853 
21810 
201519 
3 23324 
323310 
304779 
324870 
324857 
32-5088 
3.6956 
3.7911 
4.2168 
4e2157 
4.1620 
423916 
4 64036 
4.5627 
4 64145 


EQUIMOLAL MIXTURE OF CO2-CH4-C2H6 
203015 


AH 
J/MOL 


340.2171 
398.765 
467.051 
466.712 
493.601 
488.439 
559% 593 
556.614 
654.700 
777.2203 
TiS eoN 2 
802.302 
193 2962 
7862499 
737.2400 
719.171 
595.154 
443.469 
443 6433 
431.559 
278.520 
29 3905 
184.033 
172.543 


—(,H/aP) 
J/MOL-BAR 


37-490 
38.000 
38.576 
38.680 
396920 
39.665 
44.692 
44.602 
51.466 
60.571 
60.566 
65.373 
65.243 
65-370 
70.040 
HT5394 
67.717 
50.845 
SO6I Zz 
50.723 
536832 
SS ts 
22.2658 
22.667 


. o 
abbS & 


ie 


: 


gSSGett 


nw 


i 


eg 


TABLE XXXV (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 


MIXTURE OF CO2-CH4—-C2H6 
TEMPERATURE (K) 


PI 
PSIA 


267.4 
21566 
275.6 
314.4 
347.26 
368.4 
406.4 
406.4 
419.6 
419.6 
420.7 
614.6 
614.6 
613.4 
715% 
718.6 
718.6 
814.6 
814.6 
819.4 
85366 
853.6 
856.4 
987 4 
987.4 
987.4 
987.4 
987.4 
987.4 
987.4 
987.64 
10574 
TD ie-t 
1121.4 
Pi2t.4 
1216.4 
1216.4 
ielo st 
131,564 


EQUIMOLAL 


KP 
PST 


124.3 
129.0 
129.0 
144.3 
160 .0 
166.0 
167.6 
167.1 
175.1 
175.0 
158.6 
179.0 
178.1 
165.5 
160. 1 
163.9 
163.5 
148.4 
147.7 
144.6 
144.5 
143.5 
129.6 
ee 
Wee acs) 
1260 
171.5 
159.0 
158.0 
138.5 
137.8 
160.5 
138.0 
i2cet 
122.3 
TB ol) 
116.4 
127.5 
127.3 


POWER 
J/MIN 


3926092 

422.358 
422.761 

5656133 
725.265 

805.086 
936.895 
929.545 
100 6.886 
10 11.034 
881.820 
1714.355 
1695.904 
1591.088 
1950.160 
2005.081 
209 4.983 
2354.089 
2328.537 
2310 429 
25536024 
25384246 
22632890 
5111.978 
50 81.935 
5077.787 
5047.374 
44822825 
4444,522 
37362590 
3729.116 
5491.198 
4528.551 
40 20.247 
499 6.903 
3271-451 
3211.285 
2581-864 
25630999 


273215 


FLOW 


MOL/MIN 


1.0333 
1.0719 
1.0719 
1.2486 
1.4099 
12.4949 
1.6590 
1.6534 
1.6905 
1.6891 
1.6309 
203176 
223076 
2eellOL 
a sy Ah 
225502 
2 65502 
rasa ess Te) 
22/276 
Zeon 
2.1980 
Pet I39 
221090 
324635 
34661 
304778 
304752 
3.2918 
32194 
3.0717 
3.0761 
325492 
che oe ta a 
Ze pals)e) 
3.5404 
4.0363 
4.0560 
4.5573 
4.5516 


AH 
J/MOL 


379.889 
394.015 
394.391 
453 .084 
514.380 
538.037 
564.712 
5626172 
595.606 
598 .547 
540.666 
739.700 
734.900 
707.108 
770.889 
(86.221 
786.182 
859.354 
853.682 
843.988 
914.387 
908 .468 
835.663 
1475 .921 
1466. 150 
1460.053 
1452.386 
1S6ino Ls 
1355.261 
1216 .434 
t2 b2ei2ot 
15472132 
1350.136 
1130.686 
1131.756 
810.487 
(91.122 
566 2525 
5636220 


A-54 


-(AH/ AP) 


J/MOL-BAR 


44.326 
44.300 
44.342 
45.540 
46.627 
47.053 
48.869 
48.794 
49.335 
49.606 
492443 
59.935 
59.847 
61.968 
69.836 
69.574 
69.740 
83.988 
83.829 
84.654 
91.779 
91.820 
93.520 

121.974 

122.351 

122.548 

122.828 

124.223 

124.407 

127.385 

127.593 

139.808 

141.899 

133.980 

134e21g 
99 .032 
98.651 
64.445 
64.170 


“a 


EST, a! 
THs ek f 
fae. TS! 
EC. TT 
AOR, BET 
a 
OBE. céI 
TIS Ae! 
SEO. SO 
lta, se 
oP, aAG 
OFLA 


wee 


«ihds 
,eaes 


A-55 


TABLE XXXV (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
EQUIMOLAL MIXTURE OF CO2-CH4-C2H6 


TEMPERATURE (K) = 273.15 
PI AP POWER FLOW AH -—(AH/ AP) 
PSIA PSI J/MIN MOL/MIN J/MOL J/MOL-BAR 
1412.4 h25..0 1674, 845 4.7000 356.348 41.347 
1412.4 124.8 1668.690 4.26935 355 552T 4¥.503 
1613.4 119.0 851.603 4.1993 177.442 2 teaeo 
1613-% ERG. 2 827.367 421622 173.734 212500 
1613.4 116.6 825.183 4.21486 V7361t3 21.615 
1817.4 205.2 567.961 4.9660 114.369 13.800 
1817.4 119.7 567.961 4.9312 bys. 75 134955 
20 33.26 106.1 323.847 4.29231 65./86 8.992 
20 33.6 105.8 323.877 4.9055 66.022 9.050 
TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
EQUIMOLAL MIXTURE OF CO2-CH4-C2H6 
TEMPERATURE (K) = 263.15 
PI AP POWER FLOW AH —(AH/ AP) 
PSIA PSI J/MIN MOL/MIN J/MOL J/MOL-BAR 
26926 122..5 442.983 1.0631 416.653 49.331 
269.6 L228 437.882 1.0587 413.566 49.125 
348.6 154.2 789.237 1 4263 553.309 52.043 
349.6 155200 808.389 1.4426 560.359 52.198 
420 .6 185.0 1256.314 127814 705.210 Jel ell t 
420.6 184.7 1249.514 eco 701.906 esl ya heh s) 
420 .6 137.8 859.040 1.5899 540.304 56.868 
517.6 173.4 1574.284 22-0747 (586 19 63.468 
517.6 172-8 1561. 740 2 0668 135.099 63-420 
517.6 L7225 1551.990 22-0635 UscelLuD 63.2236 
620 .6 209-4 2750 847 2 6055 1055.779 PhS th care 
620.6 209.4 2750 .005 226055 1055 .456 732104 
620.6 144.0 17193275 2 02386 767.998 77.353 


(9A NH AI— | 


RAB JOM L, 


Bet . e+ 
aS 1.P+ 
£40.92 
#Of «Se 
yas .22 
Al{ ee 
868.0e 
Bde Ea 
084.4 
OES .EA 
PSL.EY 


eee tt 


TABLE XXXVI 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
MIXTURE OF CO2(.25)-CH4(.5)—-C2H6(.25) 


TEMPERATURE (K) = 363.15 
PI AP POWER FLOW AH -(AH/ AP) 
PSIA PSI J/MIN MOL/MIN J/MOL J/MOL-BAR 
Zie5 17.20 2752140 1.1996 229.342 19.339 
a1 eS 171.0 Die 209 1.1952 226.963 19.250 
39735 170 .0 268.864 1.1907 225.784 19.263 
422.5 185.1 378.732 Pe5232 248.635 19.482 
614.5 203.8 574.947 2.0673 278.103 19.791 
614.5 20 322 572.480 2.0661 277.080 1S. WUE 
614.5 202.7 569.316 2.0571 276.146 19.802 
817.5 204.3 70 8.245 2 04847 285.036 20.235 
1009.5 196.0 771.254 2.7998 275 461 20.383 
1215.5 196.3 875.964 3.1595 2774242 202484 
1215.5 196.0 874,591 3.1595 276 807 20.483 
1215.5 195.5 869.678 3.1579 275.395 20.431 
1415.5 184.1 871.225 3.3760 258.063 20.330 
1415.5 182.5 861.078 3.3607 256.218 20.362 
1617.5 187.8 966.356 3.3221 290.883 222464 
1617.5 187.5 960.948 3.7086 259.111 20.043 
1815.5 165.7 842.595 3.7504 224.665 19.665 
1815.5 164.5 829.216 3.7354 221.986 19.5%2 
2049.5 155.5 783.325 3.8998 200 .860 18.734 
2049.5 154.9 778.882 3.8933 200.056 18.731 


20 49.5 147.2 720.927 6 [oie )T/ 190.028 18.723 


(9A \HAI- 
H4B-JOMAL 


CEE. CS 
025.0! 
“88 ef 
eT .el 
T¥T.¢! 
S08, C! 
GE8405 
£8E.05 
282.05 
EH, OS 
Eres 
OFF LOS 
SAE LOS 
Soe, 6S 
FAG eS 
#38 .G4 
SsTe.e! 
tT. Sf 
Fev ,8! 
EST. BE 


HA 


e04 ats 


A-57 


TABLE XXXVI (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
MIXTURE OF CO2(.25)-CH4(.5)—-C2H6(.25) 


TEMPERATURE (K) = 333.15 

Pal AP POWER FLOW AH —(AH7 AP) 
PSIA PSI J/MIN MOL/MIN J/MOL J/MOL-BAR 
ZOO ie) 141.8 230.901 1.0221 225.904 23-106 
320.5 166.8 333.096 1.2385 268.948 23,5030 
428.5 1 Ze'O 512.868 1.6301 314.617 2340972 
608.5 19360 689.526 221132 326.290 24-406 
608.5 193.0 685.149 221094 324.793 244395 
608.5 193.0 682~ 959 221119 3236343 24.301 
S12.> 193.0 867.583 225889 335.108 25.2183 
BO 6D 192.4 863.495 2 25869 333.789 25.162 
1009.5 187.9 988.240 229498 335,045 206 O09 
1909.5 187.7 986.245 209472 334.636 256851 
1009.5 187.7 984.511 229463 334.148 25.820 
PAZ 61> 185.3 Lit2.4d43 33127 335.905 26.291 
Perks ies) 163.0 977.478 322866 297 2407 26 «463 
13.5.5 162.4 970.033 Bee oet 2952495 26.390 
1415.5 166.0 1057.050 304/798 303.760 26 2540 
1415.5 165.5 1047.206 304754 301.314 26.406 
1565.5 161.0 10 70.825 326664 292 062 26.310 
1565.5 160 .0 1056.108 3 «6538 289.040 266208 
156525 158.6 10 44.875 326398 287.069 26.252 
1715.5 159.5 10 95.808 3 68893 281.746 252620 
1715.5 159.4 1091.680 328919 280 6497 Zoe vec 
1871.5 149.5 1013.855 39836 254.2504 24.690 
1871.5 148.0 1004. 240 329631 2S% OU 24.832 
118 JAL otD 146.5 982.897 3.9413 249.383 24.689 
1871.5 140 .0 916.900 328524 238.001 242656 
203 Ue 3D 132.65 B36e207 3.9498 211.720 ERS LT Be) 
20 37.25 132.4 Bat e429 39490 212.058 23.230 
20 alee 131.8 $325 146 39426 211.065 230220 


20 37.25 DSi20 S23 e512 Shy sienhe) 209 2442 23.188 


ARI LES 


TABLE XXXVI (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
MIXTURE OF CO2(.25)-CH4( .5)-C2H6(.25) 


PI 
PSTIA 


321.4 
415.4 
421.4 
421.4 
421.4% 
616.4 
813.4 
100764 
1007.4 
1007.4 
1109.4 
1109.4 
1216.4 
1316.4 
1413.4 
1563.4 
1563.4 
1717.4 
1717.4 
1867.4 
1867.4 
1867.4 
29 39-4 
20 39.4 
29 39.4 


TEMPERATURE (K) = 


AP 
PSI 


162.3 
184.7 
187.4 
187.0 
186.6 
HOO eal 
195.5 
190.8 
190.3 
190 .0 
166.1 
165.6 
1580 
159.0 
154.1 
160.1 
159.8 
145.8 
145.2 
136.0 
ie eee 
129.4 
LES. 
117.8 
116.9 


313.15 
POWER FLOW 
J/MIN MOL /MIN 
383.445 1.2766 
568.778 1.6322 
588.981 1.6579 
586.830 1.6579 
582.791 1.6547 
Baheiii2 222109 
110 3.845 21243 
13036296 321434 
12966973 3.1401 
1288 .975 321335 
PV6i2.725 341582 
153.29 36 1524, 
1178.006 Be2g D0 
1259.920 324855 
U270.359 326156 
1429.210 329560 
1421.867 329534 
V2700 913 4.0355 
1267.520 4.0266 
1138.564 4.1200 
1190.317 4.0767 
10 56, 330 400254 
891.980 4.0829 
876.605 4.0593 
863.854 4.0342 


AH 
J/MOL 


300.342 
348 .454 
2Doe2 oe 
BO 3i67 95 
352 195 
376.2185 
405.181 
414.606 
413.027 
411.344 
3687156 
36d 60812 
3572425 
361.469 
35. Ne305 
361.272 
359.090 
314.947 
314.779 
2762346 
269.898 
262 6414 
218.462 
215.948 
214.130 


A-58 


-(AH/ AP) 


J/MOL-BAR 


260839 
27.362 
272494 
272452 
2ile3t4 
28.610 
30.059 
S) Wee ewe) 
31.479 
31.400 
322147 
32.039 
322665 
Eyer M2 
33.069 
322/28 
32.642 
31-330 
31.442 
CASTPEL ST |) 
29.2499 
29.2412 
262537 
26.588 
May yale 


aig) 
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TABLE XXXVI (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
MIXTURE DF €02(525)=€H4( 65 )=€2h6(225) 


TEMPERATURE (K) = 293615 
PI AP POWER FLOW AH -(AH/ AP) 
PSIA Pst J/MIN MOL/MIN J/MOL J/MOL-BAR 
261.4% 299.416 1.0623 281.839 S0 e702 
32744 456.306 1.3377 341.091 312114 
419.4 690.714 1.7116 403.527 305999 
419.4 719.774 1.7086 421.2263 332442 
613.4 1184.079 2 «4013 493.087 34.465 
812.4 1526.031 229312 520.607 37-510 
1910.4 1743.2 73 3 23389 522% 105 40.932 
1010.4 1734.217 323342 5206145 40.909 
1212.4 P9826541 3.7502 O25 4032 43.887 
1212.4 1971.079 367421 526.729 43.930 
1312.4 999.207 3.8391 497.302 44.799 
1312.4 1900 .309 328364 495.327 44.705 
1413.2 1947.715 4.0255 483.834 44.754 
1512.4 1892.049 4.1667 454.084 43.760 
1512.4 1881.796 4.1582 452.548 43.757 
1663.4 T55G6eel2 4.2150 369.6215 40.233 
1663.4 1547.683 4 22066 367.915 40.212 
1815.4 1450 446 464636 324.948 352489 
1815.4 1422.072 4 64433 320.048 3543853 
1815.4 1420 .996 464365 320.296 35.380 
181564 1362602 4.3851 310.677 35.208 
20 28 64 905.988 4.3598 207 803 20.01 52 
29 28.4 893.769 4.3413 205.815 2%~0L45 
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£PO, Ise 
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is er te 


THO. LR* 
Toa0OS2 
204. S82 
ALE gOR2 
£ES.A58 
eS. 454 
fac, 17s 
y Se ane 
oEA LE Ae 


Pee Bee 
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. $e. Yar 


Bae, ASE 
seG.05€ 
YTO.HLE 
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TABLE XXXVI (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
MIXTURE OF CN2(.25)-CH4( .5)—C2H6(.25) 


TEMPERATURE (K) = ZdelD 
PI AP POWER FLOW AH (7 PP) 
PSIA Bs J/MIN MOL /MIN J/MOL J/MOL-BAR 
20 Zed 29 56 367.182 1 21326 324.177 366223 
Cie ye. Er (ere 51 Dahl 9 124409 402.352 Bigewtiee 
6h) 165 2) 156.6 572.060 1.4336 399.031 36.957 
420-5 190.6 922.2269 1.8427 500 .473 38.083 
420.5 190.5 928.260 1.8456 502.944 BGrecos 
61325 179.8 1315.568 24195 543.721 43.859 
Bi) 335 S539 20 10. 733 3.0739 654.110 BeOS 
Ba 185.2 L997 25> 320672 651.189 BO e994 
IDS 395 183.2 2857.021 Bhiisre ress: U0 e)2) 61.142 
10,13 35 182.7 2843.097 326938 769.685 61.102 
1115.5 171.5 30 63.815 3.9089 783.196 66.285 
1115.5 170.8 30 48.990 328984 782.094 66-412 
B25 6 147.4 2805.575 329607 708.346 69.2699 
1215.5 147.0 2792.617 329562 TOS eat 1 69.645 
1325 61 154.2 3179.131 4.3694 727.2581 682435 
1411.5 Lys 30 28.045 4.26403 652 2546 62.603 
1411.5 150.7 30 11.289 4.6263 650.904 62.644 
1508.5 142.2 2486.3/70 421646 521-840 Ss Gy a7 2) 
1508.5 141.3 24622425 4.1488 bUSisio 2 DIZ 4 
1659.5 143.9 1976.639 Siva beh] 386.369 38.942 
UiZ:1.5 123.0 1140 .630 4 9633 229.808 262988 
1821.5 | ara) Iie tl 46 429448 LESS oh 26-940 
1821.5 119.1 1071.676 4.8740 219.871 Zo~et 15 
20 30.5 117.1 fi2Oge oie 5 O20 143.336 Ligrs3 
20 30.5 116.5 717.082 5 20204 142.833 lWieis2 
20 30.5 116.0 713.062 429986 142.649 17.835 


a 
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tr 


(AA\H Al 
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ESS «GF 
SS EeTE 
Tee .@F 
£€80.8£ 
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TABLE XXXVI (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
MIXTURE POR SCO 22 C4 25 C2H64ie25,) 


TEMPERATURE (K) = 263.15 

PI AP POWER FLOW AH =(,H/AP) 
PSIA PSI J/MIN MOL/MIN J/MOL J/MOL-BAR 

2 either 124.5 381.051 1.1260 338.381 39.420 
324.7 148 .0 580.372 1.4111 411.287 40.305 
423.7 158.3 841.294 Pantie 473.361 43.370 
423.7 1 5b 841.294 1.7811 472.332 43.330 
615.7 LTO 1477. 360 264644 599 478 5 PALS 
615.7 170 0 1476.989 2 4644 599.327 51.132 
Gieien MOS 23752540 3.1337 758.048 64.598 
815.7 171.0 2387.989 3.1370 761.214 64.564 
912.7 Toes 2602.301 3.3041 787.589 75.300 
912.7 Taleo 2590 .558 3.3017 784.595 75.361 
1009.7 161.0 3600 .626 3.7550 958.881 86.381 
1009.7 159.7 3563.855 3.7509 950.129 86.289 
1009.7 159.5 3559.091 3.7419 951.122 86.488 
1081.7 144.6 3578.233 3 68275 934.861 93.769 
1139.7 he We 4144.944 4.1463 999 .665 95.639 
1aZiet Dare 4489,967 445638 983.820 90.597 
1220a1 156.7 4455537 465505 979.126 90.625 
lation 147.3 36532621 4.8373 755.291 74.369 
1413.7 140 .0 2650-769 5.0314 526.837 54.579 
Ve isirn 13285 149 1.536 5.2578 2662560 29.178 
1818.7 12503 815.006 5.3184 153.240 120% 3% 
2021.7 Tied 468.467 5 61552 90.871 11.746 
20 2467 Lite 458.927 5.1340 89.388 11.679 
2024.7 110.0 454.335 5.1035 89.023 11.738 


7 


alt 
met 


A-62 


TABLE xXxXxvVII 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
EQUIMOLAL MIXTURE OF CO2-CH4-N2-C2H6 


TEMPERATURE (K) = 333.15 

PI AP POWER FLOW AH <A /RPA 
PSIA J/MIN MOL/MIN J/MOL J/MOL-BAR 
415 6D 429.245 125222 281.983 194.9392 
415.5 425.893 bedi 9d 280.352 19_. 922 
415.5 4176 V1 12-5191 274.597 19.915 
415.5 49 6.962 1.5004 24 ie2esd 19.898 
Dill eD 449.715 1.7275 260.319 20.061 
617.5 543.574 1.9835 274.039 20.299 
810.5 623.311 221822 205% 69 1 20.662 
810.5 622.585 21800 285.578 20-668 
1015.6 705.873 204982 282.543 20.822 
1215.6 765.084 2/631 276.890 20.797 
1414.6 798.078 209775 268.035 208579 
1414.6 789.815 2.9768 265.2315 20.555 
1414.6 780.994 229666 263%205 20-472 
1614.6 784.009 31571 248.328 20.098 
1614.6 7772186 321658 245-491 20 .048 
1614.26 770 6474 3 21484 244.717 20.030 
1813.6 693.393 323456 2072252 19.306 
1813.26 680.825 304837 195.427 195255 
1813.6 659.660 323714 195.660 19.384 
1813.6 673.077 304674 194.115 19.244 
1813.6 664.880 324542 192.481 194'252 
1813.6 648.452 304144 189.913 19.168 
202326 6322434 326306 174.194 18.059 
20 23.6 623.708 3 26086 172.835 18.099 
20 2326 587.434 325323 166.299 18.121 


(QA\HAl= 
AAG-jJGM\E 


Ss?.e! 
aTT.ef 
8ses.ef 
140.08 
eet, os 
$38, 05 
8ea.08 
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TABLE xXxxVII (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
EQUIMOLAL MIXTURE OF CO2-CH4—-N2-C2H6 


TEMPERATURE (K) = 303.15 
PI Ap POWER FLOW AH -(AH/AP) 
PSIA PSI J/MIN MOL /MIN J/MOL J/MOL-BAR 
41564 190.1 480.721 1.4711 3260776 24.931 
416.4 190.3 482.089 1.4711 327.706 24.976 
614.4 191.0 663.794 1.9433 341.578 25.938 
614.4 190.1 660 .030 129457 339.222 25.881 
81364 183.2 786.361 2oSUs 339.711 26.894 
813.4 182.2 780.798 Phe 7| 337.600 26 2874 
1013.4 189.5 986.420 Aets2T 360.960 272626 
1211.4 190.1 1131.069 3 OTT 367.501 28.038 
T2 eps 189.6 1127.485 3.0747 366.692 28.050 
1211.4 189.1 11272482 3.0747 366 691 28.124 
1413.4 180.5 1158.717 363275 348.224 27.981 
1413.4 180.2 1152.062 3.3187 347.140 27 «940 
155764 167.2 10 70 «390 3.3899 315.758 27.390 
155704 166.2 10 63.394 3.3925 313.445 27.353 
1710 4 161.8 10 20.516 NeeOUT 292.597 26.228 
1710.4 160.3 10 10.123 304877 289.617 26 «204 
1710.4 158.5 997.124 3.4851 286.104 26.180 
186704 146.8 868.417 324864 249.081 24.609 
1867.4 145.5 858.415 3 04550 2482452 24.766 
20 27 4 143.7 ed lies 325163 223.863 22.594 
202724 141.0 769.763 3 04968 220.128 226643 
20 27 4 140.8 769.763 3.5072 219 476 22.608 
20 2764 139.5 757.690 3.4968 2162675 22.527 
20 27 «4 136.7 738.706 304753 212.553 22.551 
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TABLE XXXVII (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
EQUIMOLAL MIXTURE OF CO2-CH4-N2-C2H6 


TEMPERATURE (K) = 283.15 

PI AP POWER FLOW AH —(A H/AP) 
PSIA PSI J/MIN MOL/MIN J/MOL J/MOL-BAR 
223-4 111.6 174.517 0.8110 215.176 27-964 
223.4 111.0 Ni3.992 0.8052 216.024 28.226 
27344 133.0 264.477 1.0117 261.401 28.506 
321.4% 150.7 359.380 1.2139 296.033 28.481 
425.4 188.2 606.380 1.5922 380.843 295920 
46664 200 «6 693.282 1.7002 407.761 29-482 
562.4 192.4 804.377 1.9747 407.338 30.658 
562.4 L92et 798. 748 1.9699 405.459 30.612 
71324 195.9 10 21.434 Zesoe) 434.260 322151 
85944 181.2 119 1.316 2 26018 423.289 3368S 
1012.5 Lae 3 1210.317 2.8683 421.961 35.314 
1163.5 163.1 1263. 216 3.0843 409.551 362419 
1260.5 15365 P23le207 301722 388.119 36 6672 
1413.5 161.0 1413. 769 3.5158 4022114 360224 
1558.5 146.7 1252.298 325848 349.327 S4e0D" 
1711.5 147.5 1225.01 321993 3232220 ST. 782 
1865.5 13724 10 10.529 361993 265.977 28.076 
1865.5 134.6 988,668 3.7993 260.223 28.040 
2011.5 140.9 945.761 329640 238.586 24.2559 


—— 
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TABLE XXXVII (continued) 


TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 


ea 
PSIA 


225-3 
269.3 
321.3 
323.3 
416.3 
420.4 
561.4 
615.4 
760 4 
912.4 
1005.4 
1161.4 
1261.4 
1261.4 
1261.4 
140 8.4 
1565.4 
156564 
1706.4 
1706.4 
1856.4 
1856.4 
201524 
20 1664 
20 16.4 


TEMPERATURE (K) 


AP 
PSI 


109.4 
128.0 
149.7 
150.6 
189.8 
196.0 


195.4 


179.2 
183.2 
169.5 
134.9 
144.8 
15762 
P5532 
150.5 
146.1 
145.3 
144.8 
142.0 
140 .0 
132.2 
131.4 
118.8 
119.2 
Gro, 


POWER 
J/MIN 


214.927 
3126493 
453.404 
458.225 
719669900 
798.788 
10 72.526 
19 64.553 
1410 .298 
1615.825 
1379.476 
1833. 884 
22024312 
2b626 11% 
20 58.315 
1966.599 
1700 .205 
1694. 263 
1395.804 


1373.656 


1018.697 
109 6.008 
681.711 
686, 333 
684.337 


FLOW 
MOL/MIN 


0.8664 
1.0551 
1.2808 
1.2881 
Lett58 
1.6800 
22-1017 
201778 
205827 
28861 
229538 
3 4423 
368259 
3.8081 
301364 
4.0306 
4.2212 
4.2180 
4.3458 
4.3458 
423264 
4.3310 
421915 
4 62433 
422106 


EQUIMOLAL MIXTURE OF CO2-CH4-N2-C2H6 
263.215 


AH 
J/MOL 


248.052 
296.156 
353.989 
855.129 
464.481 
475 2448 
510.301 
488.818 
546.039 
559.678 
467.014 
D326 "36 
575.647 
567.911 
550.873 
487.2915 
402.774 
401.668 
3216 Lo 
316.083 
235.2460 
232629 
i era tee he) 
161.743 
W62.524 
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(VK H7AP) 
J/MOL-BAR 


32.589 
65.057 
34.296 
34,259 
35.493 
35.182 
37.877 
29 .Das 
43.4229 
47.2890 
50.211 
53.36) 
53.111 
53.072 
53.088 
48.436 
40.204 
40.232 
32.805 
322745 
Lowe 
25-638 
19.059 
19.680 
19.858 


183 
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TABLE XXXVII (continued) 
TABULATED EXPERIMENTAL ISOTHERMAL DATA FOR 
EQUIMOLAL MIXTURE OF CO2-CH4—-N2-C2H6 
TEMPERATURE (K) = 243.215 

PI NP POWER FLOW AH =—(NHAAP) 
PSIA PSI J/MIN MOL/MIN J/MOL J/MOL-BAR 
22304 104.0 255,639 0.9013 283 2622 39 55S 
269.4 122.0 377.421 1.1046 341.653 40.616 
319.4 143.2 548.625 1.3351 410.917 41.619 
365.4 161.0 62r ks eels) 1.5612 481.252 43.353 
439.4 176.3 10 37.230 1.8588 Sores) 45.905 
439.4 121.0 623.234 Lie 938 39 Oe! 46.870 
1310.8 135.0 eM S Ne UM are 342608 355 0487 38.191 
1310.8 Ne) ere) i127.835 32622 345.727 37.617 
1413.8 127.0 1000 .310 322009 312.505 556069 
1413.8 124.2 BI 6eo29 3.1994 AOC 0iT/ 35.653 
1563.8 ibe snes Th 756.385 323628 224.924 24.770 
1714.8 P20 5226315 3 4583 151.030 17.248 
1863.8 145.0 483.065 367040 130.416 13.045 
1863.8 143.0 480.100 oof LOD 129.401 13.124 
2015.8 133.0 321.090 325498 90.452 9.863 


2015.8 131.8 317.981 35063 90.687 SRN) 


